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Abstract 

We compute exact thermodynamic geometric properties of the non-abelian quarkonium bound states 
from the consideration of one-loop strong coupling. From the general statistical principle, the intrinsic 
geometric nature of strongly coupled QCD is analyzed for the Columbic, rising and Regge rotating regimes. 
Without any approximation, we have obtained the non-linear mass effect for the Bloch-Nordsieck rotating 
strongly coupled quarkonia. For a range of physical parameters, we show in each cases that there exists 
a well-defined, non-degenerate, curved, intrinsic Riemannian manifold. As the gluons become softer and 
softer, we find in the limit of the Bloch-Nordsieck resummation that the strong coupling obtained from the 
Sudhakov form factor possesses exact local and global thermodynamic properties of the underlying mesons, 
kaons and D„ particles. 
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1 Introduction 



Quarkonia are the central objects to explore the non-perturbative nature of non-abelian gauge theories pQ. 
The fate of the heavy quarkonia in hot QCD medium describes the nature of confined-deconfinement phase 
transition (2HI] in QCD and thereby the nature of the inter-quark forces [5j-[Z]- An understanding of such 
issues is essential towards the formation of hot and dense states of quark gluon plasma matter in the heavy ion 
collision [8] and in the early universe [9]. Thus, our prime focus is to address the stability issues of vacuum 
constituted by various quarkonia states. This is in direct concern with the fact that the quark matter formation 
is an important state of art in the particle physics and its behavior is an intense topic of experimental research 
in the hot QCD media jlOj . This motivates us to investigate the nature of stability when the quarkonia possess 
a non-zero mass and spin. Following the notions of the quarkonium vacuum, we notice that an extensive study 
is required for a proper understanding of the thermodynamic stability and global statistical correlation in the 
non-abelian gauge theories. To the best of our knowledge, this is perhaps the first attempt in this direction. 

It is by now well established that the thermodynamic intrinsic geometric examinations are important in 
the configurations involving black holes in string theory [TTHT7] and M-theory [TBT - I2T] . possessing a set of rich 
stability structures [221426] . There exists a wealth of literature where thermodynamic geometry plays an im- 
portant role towards the understanding the physics of black-holes, black branes, black rings and various higher 
dimensional objects. Attentions were paid on the equilibrium perspective of black holes, and thereby one expli- 
cates the nature of the pair correlations and associated stability of the solutions. Recently, the thermodynamic 
geometry has emerged as one of the most powerful tool in the understanding covariant statistical correlations 
and the limit of the underlying microscopic-macroscopic duality issues for the 1/2- BPS black holes in string the- 
ory |26j . There have been several general notions explored in the condensed matter physics |27H34j , motivating 
us to consider the quarkonia with a set of equilibrium parameters and analyze the underlying parametric pair 
correlation functions and global correlations. We find thus that the intrinsic geometric investigation possesses 
an exact mathematical nature of the massive fluctuating quarkonia. In this concern, we provide a brief account 
of the concerned strong interactions and thermodynamic geometry in the next section. 

Given a definite covariant intrinsic geometric description of the quarkonia, we examine (i) the conditions 
of thermodynamic stability of the configuration, (ii) its parametric global correlation function in terms of a 
chosen set of fluctuating parameters. In this analysis, one can enlist the complete set of non-trivial parametric 
correlation functions of the quarkonia. The similar considerations remain valid over the black hole configura- 
tions in general relativity [551438) . attractor black holes |3T)l - fl5] and Legendre transformed chemical configura- 
tions (44][45]. Interestingly, we have explored the issue of thermodynamic geometry to the hot QCD [46] and 
offered the detailed picture of such duality relations between the thermodynamic geometry and quark-number 
susceptibilities. The issue of the quark matter stability is quite novel in its own right and calls for further 
investigations of the thermodynamic geometric structure of the phases of non-abelian gauge theories and quark 
matter production [47] . 

The thermodynamic geometric consideration is examined for the strong coupling of massless quarks which 
reveal limited effects of the QCD calculations, which have been known for the gluons since decades [481149] . 
The intrinsic geometric analysis of the linear Regge rotating quarkonium configuration [50j offers the thermo- 
dynamic stability of quark matter with a precise account of non-zero Regge mass. This consideration supports 
the formation of the topological Heiscnbcrg spin glass |51j and vector boson pair production at the hadron 
colliders [52] . Further, the thermodynamic geometric properties take an intrinsic account of the proton-proton 
and proton-antiproton scattering, hadronic total cross sections through the soft gluon summation in impact 
parameter space [53], and production of the gluons [54] . Thus, the differential geometry plays an important 
role in the study of the physics of quarkonia. 

In this paper, we examine the thermodynamic stability properties of the quarkonia system described by two 
and three parameters. We find that our analysis uniformly applies to all possible QCD particles with varying 
masses from the range of the charmonium states to the bottomonium states. This follows from the fact that the 
effective mass of the theory is kept very arbitrary; therefore the finite temperature behavior of the quarkonia 
is captured through the present analysis. From the thermodynamic geometric viewpoint of the hot QCD [46] . 
it seems quite interesting that we can reveal the near equilibrium behavior of quarkonia over an ensemble of 
hot QCD states. To do so, we introduce the thermodynamic metric tensor through the variations of the strong 
coupling, and focus our attention on the thermodynamic properties when gluons become softer and softer [53] . 
Our consideration reveals the thermodynamic geometric nature of very soft gluons and thereby describes the 
stability properties of non-abelian quarkonium bound states with and without the mass. 

As explained in the section 2, the problem thus posed starts from the very definition of Sudhakov form 
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factor [55], and thus the intrinsic geometric understanding of the QCD coupling in the limit of the vanishing 
transverse momentum energies. Notice that our analysis considers all types of QCD states, as we do not 
truncate the theory but work with the full Bessel function underlying in the Sudhakov form factor. We initiate 
our analysis where the computation requires only the first few excited states in which one may think that only 
the charmonium states play the dominant role in the quark matter formation. The analysis is continued for 
the mass-less Regge rotating quarkonia, and finally all possible mass states are considered through the strong 
coupling. Our study is two folded, (i) the investigations on the zero temperature properties of quarkonia, 
and (ii) extensions to the finite temperature, as well as very high temperature properties of the quarkonium 
vacuum fluctuations in the non-abelian gauge theories. This shows the universality of our consideration of the 
thermodynamic geometry, that applies for all possible mass excitations of the quarkonium vacuum, and thus 
the stability of quark matter formations. 

The effect of finite mass is explored for the Regge mass and arbitrary mass quarkonia with the one-loop 
Bloch-Nordsieck resummation |56j . Following the Polyakov arguments |50j . our analysis shows that the rising 
potential vacuum can be thought of as a collection of decaying and weakly interacting particles. However, the 
Regge rotating configuration offers almost the same stability structure and is also nearly weakly interacting. 
Subsequently, we show that the global thermodynamic correlation of the Regge rotating states remains almost 
the same as that of the massless states, up to the sign of the rotation term. With the large quark mass, the strong 
coupling acquires a large twist and thus thermodynamically the configuration becomes highly unstable. For 
experimental reasons, the graphical characterizations are offered for large QCD cut-off, where the quarkonium 
vacuum attains the Regge rotating regime. Such an investigation clarifies the thermodynamic stability properties 
of the quark matter formation. With this motivation, we shall focus our attention on the (i) quark masses and 
their dependence on the decay constants, (ii) soft gluons and meson decays, (iii) very soft gluons and (iv) strong 
coupling hadronic form factor. 

The paper is organized as follows. In section 2, we present a brief account of the strong interactions present 
in the quarkonium formation and its thermodynamic geometry. In section 3, we consider the non-rotating 
massless quarks and examine their thermodynamic geometric stabilities. Subsequently, we extend the analysis 
for the non-zero angular momentum. In section 4, we explore generic cases in the framework of the Sudhakov 
form factor in the limit of Bloch-Nordsieck resummation for the massive quarkonia, with and without rotations. 
In section 5, we discuss some concluding issues towards the thermodynamic geometric picture of quark matter 
formation and the stability of underlying QCD particles. 

2 Phases of the Strong Interaction 

In this section, we shall set-up the formulation of problem and outline the notion of the intrinsic thermodynamic 
geometry. The non-perturbative gluonic effects gives the relation of effective QCD theories with the index of the 
potential, momentum scale, mass and angular momentum, if any, and thus the consideration for the analysis. 
Most efficiently, our intrinsic geometric model is designed to provide the critical values for the parameters of 
the gluonic effects. It is worth mentioning the importance of the present method that the remaining values of 
the parameters would be an stable mode for the considered type of quarkonia. 

2.1 Quarkonium Bound States 

In the physics of quarks, the non-perturbative gluonic effects are addressed via the considerations of Bloch and 
Nordsieck [56). Some quantities of known interest are intrinsic transverse momentum of Drell-Yan pairs [57], 
impact parameter distribution of partons and rise of total cross section [58j , form factors and decay constants }59| 
of the particles. A gluon is said to be soft, if the transverse momentum k± is kept small. In this limit, the 
infrared limit of the charmonium is always unobservable. This follows from the fact that a momentum integration 
is required to perform, if one wants to keep a vanishingly small transverse momentum, i.e., k± — > for the 
gluons. In the case of abelian theory, the results follow directly from the properties of the underlying Poisson 
distribution and their numerical counterparts. 

In the case of non-abelian theory, the consideration of Sudhakov form factor turns out to be an efficient tool 
for the soft gluons, which is described as below. Let us recall that the QCD coupling a s {k±) never lies near 
the limit k± — > 0, and thus the QCD effects are limited for the bound state thus formed after the resummation. 
This initiates the requirement for the integration over the k± and thus the necessity of Sudhakov form factor. 
Let us illustrate the computations for decays processes associated with the pions, viz., ir° — » 77 and ir — > jiv. 



2 



During the process, when a virtual photon hits the quarks, soft gluons are emitted with k± — > and a 
non-zero total momentum k 2 ^ 0, satisfying the standard transverse momentum Poisson distribution d P(k±). 
The associated Poisson summed vertex and the transverse limiting momenta are respectively given by 

d 2 P(k ± ) 



r ™ " dk\ 



n(0) 



d 2 P(k ± ) 



(i) 

d 2 t exp(-h(b)), 

(2) 

where the exponent defined as 

h(b) = [ d 3 n k (l - exp(-it .kl)) (3) 



is known the Sudhakov form factor. Recall that one is required to take the vanishing transverse momentum 
limit for a decay, which is done systematically with the help of the Sudhakov form factor, see for a review [60] . 
We now illustrate the case of masslcss quarks. For the quark mass m q = 0, in turns out that the Sudhakov 
form factor reduces to the following integral 



dki dk± /1 , >r -> 



h(b)= I (4) 



where k± € (0, m P / 2 ) is due to a physical reason that on average each soft gluon can take as much as half of the 
initial center of mass energy. Herewith, the consideration of Sudhakov form factor comes for a model of limiting 
transverse momentum QCD coupling ot a (kx.)\k±-yo- Inspired by the behavior of relativistic wave functions near 
the origin for a QCD potential [6lJ|62] f° r the Rechardson potential for the quark bound states, one arrives at 
the transverse momentum dependent strong QCD coupling 

tu2\_ 12n P (x\ 

as( ^-(33-27V / )ln(l +P ( ] ^-) 2 p) (5) 

in the limit of one loop gluon exchange potential. Towards the determination of the index p, an interesting 
argument follows from the consideration of Polyakov |50| which we shall explore further from the perspective 
of thermodynamic geometry in the subsequent discussion of the present paper. Before doing so, let us consider 
the joint effects of the (i) confinement and (ii) rotation, and make a platform to describe the properties of the 
thermodynamic geometry. In this case, the configuration is described by Regge trajectories with the leading 
order effective potential 

V(r,J)= J(J 9 +1) +Cr 2 *-\ (6) 
Such an effective potential is inspired from the limiting QCD strong coupling 

In (l+p (J*— )*) 

with the matching of the pre-factor as b := (33 — 2Nf)/12iv. In the momentum space, the net effective potential 
offers the right quarkonium bound states, after taking account of the one-loop exchange terms. In the sense 
of one- loop quantum effects, we propose that the "quantum" nature of quark matters follows directly from the 
consideration of the thermodynamic intrinsic geometric potential of Richardson types. The index p defines the 
nature of the potential |50) . i.e. whether the quarkonium lies in the Coulomb phase (p = 0) or in the rising 
phase (p > 1/2). Notice further that the Regge behavior may also be determined from interpolating values of 
the index p, which is the matter of the subsequent sections. For the future purpose, it is worth mentioning that 
the minimization of the effective potential containing the confining and rotation effects yields the value p = 5/6 
for the index, which corresponds to the linear Regge regime. We shall exploit these facts, while dealing with 
the thermodynamic geometry of rotating (massive) quarkonia. 
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2.2 Thermodynamic Geometry 



In this subsection, we offer a brief account of the essential features of thermodynamic geometries from the 
perspective of the application to the strongly coupled QCD. We shall focus our attention on the thermodynamic 
geometric nature of quarkonium bound states with a finite number of parameters carried by the effective field 
theories. 

Let us consider the framework of the intrinsic Ricmannian geometry whose covariant metric tensor is defined 
as the Hessian matrix of the QCD coupling, with respect to a finite number of arbitrary parameters carried 
by the soft gluons and quarks. This consideration yields the space spanned by the n parameters of the strong 
QCD coupling a s , which in the present consideration, exhibits a n-dimcnsional intrinsic Ricmannian manifold 
M n . Following the notion of [24Tf32] for the thermodynamic geometry, the components of the covariant metric 
tensor are defined as 

_ d 2 a s {x) 

9l] •" dx>dat ' [ ' 

where the vector x € M n . In the quarkonium effective configuration, there are only a few in physical 
parameters which make the analysis fairly simple. As mentioned in the foregoing section, the variables of the 
interest in the present study of quarkonium are the momentum scale parameter, Q 2 := q, the mass M and the 
angular momentum J, if any. Thus, let us first illustrate the consideration of thermodynamic geometry for the 
non-rotating configurations, viz., J = 0. In this case, we can take the QCD coupling a s (q,p) as a function of q,p 
and thereby may explore the local and global thermodynamic properties towards the stability of the quarkonia 
on gp-surface. For a given QCD coupling a s (q,p) := A, an equilibrium quarkonium configuration is achieved at 
the points, where the first derivatives 



OA^p) = d_A^p) =Q 
dq dp 

vanish identically, implying the existence of the equilibrium data {qo,po} on an intrinsic Ricmannian surface 
(M2 (R) , g) . In order to verify whether the quarkonium configuration is minimally coupled at {qo,Po}, we may 
follow the set-up of thermodynamic geometry and define the components of the thermodynamic metric tensor 
as 

d 2 A d 2 A _ d 2 A 

9qq ~ w 9qv ~~oWv 9pp ~W 

(10) 

In the present case of the (M2(R),g), it follows that the determinant of the metric tensor is 

\\g(q,p)\\ = A qq A pp -A 2 qp . (11) 

Explicitly, we can calculate T^k, RijkU Rij an d R for the above two dimensional thermodynamic geometry 
(Mi , g) of the non- rotating quarkonia and may easily see that the scalar curvature is given by 

R{QiP) = ~ 2\\g(q ~pjv&(~ / ^ pp '^ qqq '^ qpp ~ ^ii^qpp ~ ^-pp-^-qqp 

~\~AqpSqqpAqpp -\- AqqAqqpAppp -(- AqpAqqqAppp}. (1^) 

As a global intrinsic geometric invariant, the scalar curvature accompanies information of the correlation 
volume of underlying quark fluctuations. The scalar curvature further explicates the nature of long range global 
correlations and phase transitions, if any, deriving from a given phase. In this sense, we anticipate that the 
set of particles corresponding to the specific decay, are statistically interacting, if the underlying quarkonium 
configuration has a non-zero thermodynamic scalar curvature. Incrementally, we notice that the configurations 
under present consideration are allowed to be effectively attractive or repulsive, and weakly interacting in 
general. For the two dimensional thermodynamic geometry |24j defined as the intrinsic Riemannian surface 
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(Al2{R),g), the relation of the thermodynamic scalar curvature to the thermodynamic curvature tensor is given 
as 

2 

R{q,P) = J^Rqpqp- (13) 

The intrinsic geometric analysis thus provides a set of physical indications encoded in the geometric quan- 
tities, e.g., the scalar curvature and possible geometrically non-trivial invariants. For a given quarkonium 
configuration, the underlying analysis would involve an ensemble or subensemble equilibrium configuration 
forming the statistical basis for the Gaussian distribution of the particle, which is ensured for any distribution 
in the late time limit. With this brief introduction to the thermodynamic geometry, we shall now proceed 
to systematically analyze the underlying stability structures for the parametric fluctuations of the massless 
non-rotating quarkonia. 



3 Massless Quarkonia 

In the present section, we shall examine the nature of the QCD coupling a s (q,p) from the thermodynamic 
stability properties. Firstly, we analyze the massless quarkonium for non-rotating configurations and then take 
take an account of the rotation in the next subsection. In the subsequent analysis, we shall adopt the notation 



3.1 Non-rotating Quarkonia 

As stated earlier, the thermodynamic metric in the quarkonium parameter space is given by the Hessian matrix 
of the strong coupling with respect to the intensive variables, which in this case are the two distinct parameters 
{q,p} carried by the quark matters. Considering Eqn.([7]), one obtains the following expression: 



A ^ P) ■■= bWTpWwr (14) 

for the strong QCD coupling. To compute the thermodynamic metric tensor in the parameter space, we 
employ the Eqn. (|10p . which leads to the following expression for the components of the metric tensor 



9 11 ~ u 2 O ' 9qp — , Q i 9pp — i o ■ l 10 7 

bq 2 rg bq r Q b r Q 

In this framework, we observe that the geometric nature of parametric pair correlations offers the notion of 
fluctuating quarkonia. In order to simplify the subsequent notations, let us define the logarithmic factor as 

l(p):=Hl+p(q/Lf). (16) 

In this case, we find that the factors in the numerator of the local pair correlation functions arc expressed 
as follows 



2p 2 (q/L) 2 P + l{p){{q/Lf - p{q/Lf + (q/L) 2 Pp), 
2(q/L) 2 Pp + l(p)(2(q/L) 2 Pp 2 \n(q/L) - ^q/Lf 
-2(q/L)^p-(q/Lr P \n(q/L)), 
2(q/L) 2 P( P + 2p 2 ln(q/L) + p 3 Hq/Lf) 
-l(j>){{q/LY + (q/LfPp + Aiq/LfpHq/L) 

+2(g/L) 2 V ln(q/L) + (q/Lfp 2 \n(q/L) 2 ). (17) 

We further notice an interesting conclusion for the denominator of the local pair correlation functions, and 
we find for the massless non-rotating configuration that the denominators {r^| i = 1, 2, 3} of all the local pair 
correlation functions take uniform value l(p) 3 cxp (2l(p)). 

Thus, the fluctuating quarkonium configuration may be easily analyzed in terms of the parameters of the 
underlying effective theory. Moreover, it is evident that the principle components of the metric tensor, signifying 



,Q 



'•22 



5 




self pair correlations remain positive definite functions. In a given QCD phase, this happens when the parameters 
{q,p} are confines in the domain 

{V := (q,p) e M 2 | > 0, n| > 0}. (18) 

Over this domain of {q,p}, it is worth mentioning that the massless non- rotating quarkonia is well-behaved and 
locally stable. The global stability is offered by computing the determinant of the metric tensor, and requiring 
it to be positive definite. Following the Eqn. fTTj) . we find further for the generic value of the parameters that 
the Gaussian fluctuations form a stable set of correlations over {q,p}, if the determinant of the metric tensor 

"■ 9 " = &VZ(p) 5 exp(3;( : p)) r ^ (19) 

remains a positive function on the intrinsic gp-surface (M2(R), g). Explicitly, we obtain that the numerator of 
the determinant of the metric tensor can be expressed as 

nf(q,p) : = 2(q/L) 3p (p + 3p 2 + 2p 2 ln{q/L) + 2p 3 ln(q/L) + p 3 ln(q/L) 2 ) 

-l(p)(2(q/L) 2p +p(q/L) 3p + 2p 2 (q/L) 3p ln(q/L) + ip 2 (q/L) 3 ? 
+4p(q/L) 2p ln(q/L)+p 2 (q/L) 2p ln(q/L) 2 + 2p 2 (q/ L) 2p ln(q/ L) 
+7p(q/L) 2p ). (20) 

The behavior of the determinant of the metric tensor shows that such massless non-rotating quarkonium 
becomes unstable for opposite values of q and p for positive n^(q,p). For generic q and p, the nature of the 
determinant of the metric tensor is depicted in the Eqn. (|19p . showing that the quarkonia become unstable in 
the vanishing limit of q or p. It is worth mentioning for a common sigh of the index parameter p and 6, that the 
non-rotating massless quarkonia is stable in the regions of the (M2(R) , g) , where the function n^(q,p) picks up 
the positive sign. 

Let us explore the functional behavior of the associated scalar curvature, in order to examine the mixing 
effects of the parameters forming the intrinsic surface. In this case, there is only one non-trivial component of 
the Riemann Christoffel tensor R qpqp . As per the definition of the Eqn. (|12[ ), our computation shows that the 
scalar curvature reduce to the following specific form 

^P) - ^rL( n n )Q +^ )Q m + n^ Q l( P r + n^l( P n (21) 

2p (Tig ) 
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Figure 2: The determinant of the metric tensor plotted as a function of q, describing the fluctuations in massless 
non-rotating quarkonia near the Coulomb regime. 
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Figure 3: The determinant of the metric tensor plotted as a function of q, describing the fluctuations in massless 
non-rotating quarkonia in the raising regime. 
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The factors of the numerator of the scalar curvature take the following expressions 



f Q = \n{q/Lf - 16(g/i) 2 V In(g/L) - 16(q/L) 2p p 7 \n(q/L) 3 

~56(q/L) 2p p 6 ln(q/L) 2 - 48(q/L) 2p p 6 Hq/L) - 4(q/L) 2p p 7 ln{q/L) 4 



-16(q/L) 2p p 7 hi{q/L) 2 - I6{q/L) 2p p e ln{q/L) 3 -64{q/L) 2p p 5 \n(q/L) 
-24(q/L) 2p p 5 ln(q/L) 2 - 36(q/L) 2p p 5 - A{q/L) 2p p 3 , 

?4 1)Q = 28p 6 {q/L) 2p \n{q/L) 2 + Up 4 (q/L) 2p + 20p 5 (q/L) 2p \n(q/L) 2 
+50p 3 (q/L) p + 40p 3 (q/L) p \n(q/L) + U6p 4 ln(q/L)(q/L) p 
+56p 4 (q/L) p \n{q/L) 2 + 8p 6 (q/L) 2p Hq/L) 3 + Ap 3 (q/L) 2p 
+ I0p 2 (q/L) p + A8(q/L) 2p p 5 + 102(g/i)V + I12p 5 {q/L) p \n{q/L) 2 
+32(g/L) 2 V ln(q/L) + 2A(q/L) p p 6 Hq/L) 2 + 32p 5 (q/L) p \n(q/L) 3 
+104(«/L)V Hq/L) + 24(q/L) p p 6 \n(q/L) 3 + 6(q/L) p p 6 Hq/L) 4 
+ I6p 4 (q/L) 2p H<l/L) + 50p 5 Hq/L){q/L) 2p , 

nf Q = -6p 3 (q/L) 2p - 8p 5 (q/L) 2p Hq/L) - 6p 4 ln( 9 /L)( g /L) 2 f 
+6p 4 {q/L) 2p + I0p 3 (q/L) p - Up 2 {q/L) p - 20p 2 - 12p 

-30 P 2 Hq/l) + 8(q/L) p P e h<i/l) 2 i(p) 2 - 20 P 3 ( q /L) p Hq/L) 
-4 P 5 (q/L) 2p Hq/L) 2 i(p) 2 - 90 P 3 Hq/L) - 36 P 3 Hq/L) 2 

-62(q/L) p p 4 l(p) 2 - 16(g/L) 2 V - 16p 4 Hq/L) 3 

-8 P 5 Hq/L) 2 - 56p 4 Hq/L) - 56 P 4 Hq/ L ) 2 

~8p b Hq/Lf - 2p 5 \n(q/L) 4 - 8(q/L) p p 5 Hq/L) 
+8{q/L) p p 6 Hq/Lf + 2(q/L) p p 6 Hq/Lf - 82p 3 
-l6p 4 (q/L) P Hq/L) 2 - Mp 4 Hq/L){q/L) p , 
n {3)Q = -lhp 3 Hq/L) + 2p 3 {q/L) 2p - 4p 4 (q / L) 2p - 8p 5 ln(q / L) 2 (q / L) p 
-9p 4 ln(q/L)(q/L) p - 2p 4 ln( g /L) 2 (g/Lf - 2p 5 \n(q/ L) 3 (q/ L) p 
+p 3 Hq/L)(q/L) p - 8p 5 Hq/L){q/L) p - 20p 3 (q/L) p + 9p 2 (q/L) p 
~hp 4 Hq/Lf - V Hq/Lf - lQp 4 Hq/L) - 19/ln( 9 /L) 2 
-4p 5 ln(g/L) 3 - p 5 + 16(g/L)"V - 6(q/L)- p p + 5(q/L)- p p 2 Hq/L) 2 

+(q/L)- p P 3 Hq/Lf + io(q/L)- p P 2 Hq/L) + 3( 9 /L)-y Hq/Lf 

+6(q/L)- p pHq/L) + Mq/L)- p p 3 Hq/L) + 6(q/L)- p - 28p 2 

-24(q/L) 2p p 4 + 7p 3 - 3p 3 ln(q/Lf + 3p 2 Hq/L) + I2p. (22) 

In the case when the Ricci scalar curvature R(q,p) vanishes, the underlying quarkonium system is found to 
be in equilibrium. Such a state of the configuration can arise with {nW Q = 0, i = 0, 1, 2, 3}, if the other factors 
of the scalar curvature remain non-zero. In the other case, when the R(q,p) diverges, the configuration goes 
over a transition. Such an extreme behavior of the quarkonia is expected to happen, when either the index p 
or the the numerator of the determinant of the metric tensor vanish. Geometrically, the intrinsic gp-surface 
becomes the flat Euclidean plane in the first case, while it gets infinitely curved in the second one. 

For the choice of QCD parameters L = 1000 and 6=1, the Fig.([T]) shows the plot of the determinant 
of the metric tensor. This plots explicates the nature of the stability of non-rotating configurations in the 
limit of massless quarkonia. The corresponding plot for the scalar curvature is depicted in the Fig. ((4]). This 
plot shows the global nature of the non-rotating configuration in the limit of massless quarkonia, under the 
effects of Gaussian fluctuations of the index p and the parameter q. Notice, for the generic massless non- 
rotating quarkonium configurations, that the limiting scalar curvature interestingly simplifies to the slice shape 
of Fig. ([T]). Further, it is worth mentioning that the corresponding peaks of the determinant of the metric tensor, 
as depicted in the Fig. (jT} indicate the graphical nature of the instability present in the general non-rotating 
massless quarkonia. Physically, the presence of peaks in the determinant of the metric tensor shows a non-trivial 
interaction in the system. 

For the regime of the Coulambic potential, the respective surface plots of the determinant of the metric tensor 
and scalar curvature are respectively shown in the Figs.© and ([5]) for p = 0.1. We observe that the stability 
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Figure 4: The curvature scalar plotted as a function of q, and the index parameter p, describing the fluctuations 
in massless non-rotating quarkonia. 
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Figure 5: The curvature scalar plotted as a function of q, describing the fluctuations in massless non-rotating 
quarkonia near the Coulambic regime. 



10 




Figure 6: The curvature scalar plotted as a function of q, describing the fluctuations in massless non-rotating 
quarkonia in the rising regime. 
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of massless non-rotating quarkonia exists in certain bands. A closer view shows that the limiting Coulambic 
quarkonia are decaying and interacting particles. For the regime of the rising potential, e.g., p — 0.51, the 
respective surface plots of the determinant of the metric tensor and scalar curvature are respectively shown in 
the Figs.(|3]) and (|6]), which shows that the limiting rising potential quarkonia arc stable and non- interacting 
particles. This follows from the fact that the determinant instability is present only for a specific set q. 

3.2 Regge Rotating Quarkonia 

In the present subsection, we analyze the nature of massless rotating quarkonia generated by the parameters 
q,p and angular momentum J. Following Eqn.©, we find that the modified strong QCD coupling takes the 
form 



A(q,p,J) = l - P +lj( J+ i). (23) 

£>ln(l +p{q/L)P) bi 

To focus on the intrinsic geometry of the present case, we chose the parameters q,p, J as the variables for 
the QCD coupling. As in the previous subsection, we may again exploit the definition of the Hessian function 
Hess(A(q,p, J)) of the QCD coupling. Considering the analysis of the Regge model, we find that the components 
of the metric tensor in the framework has the same characterizations for the g qq , g qp and g pp as obtained for 
the corresponding rotating massless quarkonia. The component of the metric tensor for the rotation parameter 
turns out to be gjj = While, the remaining components of the metric tensor, involving J and either the q 
or p, vanish identically. 

It follows that the pure pair correlations {g qq , g pp ,jj } between the parameters {q,p, J} remain positive same 
as in the case of the non-rotating quarkonia. Further, our computation demonstrates the over-all nature of the 
parametric fluctuations. In fact, we find that the determinant of the metric tensor reduces to the following 
simple expression 



11511 b 2 b iq 2 l(p) 5 exp(3l(p)) U97 (24) 

where the n^(q,p) remains the same as for the non- rotating case. It is worth mentioning that the Regge 
rotating massless quarkonia is well-behaved, as long as the corresponding rotating massless quarkonia remains 
so. Over the domain of the parameters {q,p, J}, we thus observe that the Gaussian fluctuations have the same 
set of thermodynamic metric stability structures, as long as b\ > 0. The observation of the metric structure 
shows that the Regge rotating massless quarkonia is stable on the gp-surface if the index parameter p and the 
function n^(q,p), appearing in the numerator of the determinant of the metric tensor have the same sign. 

Furthermore, we may easily analyze the underlying important conclusions for the specific considerations 
of Regge rotating massless quarkonia. As in the case of the non-rotating massless case, the global nature of 
the scalar curvature and associated phase transitions may thus be determined over the range of parameters 
describing the Regge rotating quarkonia of interest. Since the previous non-rotating configuration is trivially 
embedded in the Regge rotating configuration, we thus find that the thermodynamic scalar curvature remains 
exactly the same as for the non-rotating massless quarkonia. The equality of the two scalar curvatures shows 
that the Regge rotation keeps the same global thermodynamic stability structures as the non-rotating massless 
counterpart. In the Regge trajectory model, we thus find an interesting conclusion that all possible local and 
global thermodynamic stability behavior of the massless quarkonia remains the same up to the sign of b\ , as if 
there were no rotation in the underlying configuration. 

The fact that the efficiency of the rotation induces a mass to the quarkonium is analyzed by considering the 
Bloch-Nordsicck rcsummation of the angular phases. As per the analysis of the next subsections, our method 
offers a non-linear generic characterization for the quarkonium configurations, which we consider to be neither 
purely massless nor purely rotating, but a strongly coupled QCD. 



4 Massive Quarkonia 

In this section, we extend the intrinsic geometric analysis for realistic quarkonia and compute the associated local 
and global thermodynamic quantities in the subsequent subsections. Following the notations of the previous 
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section, a given generic quarkonia can analytically be easily analyzed about some local equilibrium, if we can 
fix one of the parameter. The logic simply follows from the fact that the underlying configurations reduce to 
an intrinsic surface, and thus are easy to compute. Such an assumption is allowed for the specific phases of 
the QCD, whether one considers the Coulomb phase or the rising phase. Our analysis certainly does not stop 
here, but we indeed extend it for the general quarkonia with {q,p, 7} all fluctuating. Let us firstly illustrate 
the cases for the two parameter configurations with either {q, 7} or {q, m} fluctuating and then systematically 
for the generic quarkonia. 



4.1 Quarkonia in QJ-plane 

In this subsection, we shall use the essential features of thermodynamic geometry, thus developed, to describe 
the quarkonia using modified strong QCD coupling, with an increasing number of parameters. Let us focus our 
attention on the geometric nature of the local and global correlations in the neighbourhood of small fluctua- 
tions, in the chosen quarkonium configurations. As per the consideration of massive quarkonia, the associated 
rcsummcd strong QCD coupling is given by 



"<«• J » = \ wTkw) H 7r& {i - ma ^ ] - (25) 

where J(V, x) is the Bessel function of the first kind of the order v. As stated earlier, the thermodynamic 
metric in the parameter space is given by the Hessian matrix Hess(A(q, J)) of the strong QCD coupling with 
respect to the variables defining the thermodynamic manifolds. In order to simplify the subsequent expressions, 
let us define the logarithmic factor of the concerned Bloch-Nordsieck rotation as 



fM :=ln( ^ + ^Zj ). (26) 
V J - V 7 - q 

In this framework, it turns out that the fluctuation nature of parametric pair correlations may be easily 
divulged in terms of the momentum transfer and rotation parameter of the underlying quarkonia. Following 
Eqn. pop . wc find, under the Gaussian fluctuations of {q, 7}, that the components of the metric tensor are 



a - E (n {0)J + n {1]J l(v)+n i2)J l(v) 2 ) 

9qq ~ 46/(p) 3 exp(2Z(p))g5/2(j_ 5)3/2^11 + n n + )■> 

a = P (n {a)J + n {1)J l(v)) 

yqP 2bl(p) 2 exp(l(p))q 3 / 2 VJ(J - q) 3 ' 2 12 12 KP,h 

(h, = - 2 b l( P )J^(J-q)V^ J - q){1 - J{ °^- (2?) 

As a result, wc find without any approximation that the factors in the numerator of the local pair correlation 
associated with the gg-components are expressed as 



= 8p\f((q/LyP) 2 ^/j^(q 3/2 -q 3/2 J(0,a^)^ y/qJ + y/qJ(0,ay/q)J), 
n[\ )J = -Ap 2 f^7~q{q/Lf{J{l 1 a^)aq 2 + q 3 ' 2 ~q^ 3 ' 2 J{Q,a^) + J{l,a^)aqJ 
-sfqJ + VqJ(0, <h/q)J) + 8p 2 (q/LnVjq 3 / 2 + J 3/2 ^J(0, ajq) 
- VJq 3/2 7(0, a^q) - J 3 / 2 ^) + 4p 3 f^J^(q/L)P(^J - JqJ% a^)J 
+q 3 / 2 J(0, a^q) - q 3 ' 2 ) + Ap 3 f ^J~q({q/ Lf ) 2 {q 3 ' 2 + y/qj(0, a^q)J 
-q 3/2 J{0, a^fq) + J(l,a^/q)aqJ - JqJ - J{l,a^/q)aq 2 ) 
+8p 3 ((q/Lr) 2 (Vjq 3/2 - VJq S/2 J(0,a^q) + J 3 ' 2 \/qJ(Q,aJq) - J 3/2 ^), 
n [2 } J = 2VJ{6q 3/2 p(q/L) p + 3q 3/2 - AJ 3j2 yfq - 3q 3/2 J(0, Oy/q) - 2J(l,a^)aq 2 ) 
+.fV J - l( J (°> ay/q~)a 2 q 5/2 - 7(0, a^/q)a 2 q 3/2 J + 27(1, a^)aqj 
-27(1, a^)aq 2 ) + 47 3 / 2 (7(l, a^q)aq + ^7(0, a^q)) 
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Figure 7: The gg-component of the metric tensor plotted as a function of q and the angular momentum J, 
describing the heat capacity of the charge in massive rotating quarkonia. 




+2p(q/L)> ) (fa 2 q 5/2 ^J^J(0 1 a^) - fa 2 q 3 l 2 ^J~qJ{0,a^)J 
-2/J(l, ay /q)aq 2 ^J -q + 2fJ(l,a,y/q)aqy/j - qJ + 4 J 3/2 J(l, a^/q)aq 
-4 J 3/2 ^q + U 3/2 ^J(Q, ay/q) - AVjJ(l,a^)aq 2 - 6V7g 3/2 J(0, o^g)) 
V((?/i) P ) 2 (/a 2 '? 5/2 V / ^9^(0,«V9) - fa 2 q 3 l 2 ^J~qJ{Q,a^)J 
+2/J(l, ay/q)aqy/ J - qJ - 2/ J(l, a^q)aq 2 yj J - q - 4</jJ(l, a^/q)aq 2 
+4 J 3 / 2 J(l, a^)ag + 4J 3 / 2 ^J(0, aVg) - 4J 3/2 ^ + 6 W /2 

-6\/jg 3/2 J(0,aV9))- (28) 
While, the factors in the numerator of the associated gj-component are given by 

nf 2 )J = 2 P 2 {q/L)P{q 3 / 2 - J 3 / 2 + J(0, o>/j) - J{0, aVJ)), 
n$ J = q 3/2 + aqJJ{l,aVj)-aq 2 J{l,aVj)-q 3/2 J(0 1 a\ r J) 

+ (q/L) p p(q 3/2 + aqJJ(l,aVJ) - aq 2 J{l, a^fj) - q 3/2 J(0, aV7)). (29) 

Herewith, we see that the geometric nature of parametric pair correlations turns out to be remarkably 
interesting. The fluctuating quarkonia may be easily described in terms of the q and J. For the configurations 
with the same sign of the index p and constant 6, it is evident that the principle components of the metric tensor, 
signifying self pair correlations, are positive definite functions in a non-trivial range q. The local stability requires 
that (i) qq- fluctuations satisfy the constraint 

ng )J + n[\ )J l{p) + n [2 } J l{p) 2 < (30) 

and (ii) J J- fluctuations be constrained to the following limiting values of the Bessel function 

J(Q,a^/q) < 1, q > 2 J, 

> 1, q<2J. (31) 

A straightforward computation demonstrates the over-all nature of the parametric fluctuations. In this case, 
we find that the determinant of the metric tensor reduces to the following expression 

1 {n^ + n^ J l{p) + nf J l{pf), (32) 



tf 8& 2 Z(p) 4 exp (2l(p))q 5 / 2 J 3 / 2 (J - qfl 
where the coefficients {rig 1 * >J , n^ J } appearing in the determinant of the metric tensor factorize as 



4 1)J = V( g /L)"(4 12 > J + n£V J p(q/Ly), 

n (2)j = n (20)J + 2 nfV J p{ q /Ly + nfV J {{q/Ly) 2 p 2 . (33) 
Explicitly, we find that all factors of the numerator of the determinant can be presented as 



nf J = %p\(q/Ly) 2 (^fjJq-^J~q + ^J^^T~qJ{Q,a^) 2 
-fq 3/2 J(0, a^f - fq 3 ' 2 + 2Jf^ + 2fq 3 ' 2 J{Q, a^/q) 
-2^J^^I~qJ(Q, a^q) + 2Jf^J(0, a^) 2 - iJf^/qJ(0, a^/q)), 
< 2)J - (Ml, a^)aq 2 - fq 3 ' 2 - / J(l, a^)aq 2 J(0, a^q) - .fq 3 ' 2 J{Q, a^) 2 
+2fq 3 / 2 J{0, a^j) + 2y/jqJ(l,ay/q)ay/j - qj(0, a^q) 
+2 Jfy/q + 2JfJ{l,a^q)aqJ(0, a^q) + 2 Jfy/qJlfi, a^q) 2 
-2 Jf J(l, a^/q)aq - 2y/jqJ(l, a^/q)ay/j -q + AVJ \fq^J 1 - q 
+±y/ly/qy/j-qJ(Q, a^) 2 - 4J/^J(0, a^q) - 8y/Jy/qy/J - qJ(0, oVa)), 
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4 13)J = /(? 3/2 + 9' V2 «V?) 2 - 2 J^J(0, a^) 2 - 2^/? 

-2q 3 / 2 J(0, a^g) + 4 J^J(0, aV«)) + ((<?/£)")(/ a^W 2 
-fqV 2 J(0, a^) 2 - fJ(l, a^)aq 2 J(0, ay/q) - jq 3 ' 2 
+2/ JV? - 2y/jqJ(l,ay/q)ay/j-q + 2/g 3/2 J(0, a^) 
+2fJJ(l,ay/q)aqJ(0, a^) - 2/ JJ(1, a ^/g)ag + 2fJy/qJ{0, a^g) 2 
+2\/jgJ(l, ay/q)ay/j - qJ(Q, ay/q) + iVly/qy/ J - qJ(0, ay/q) 2 
+A\/ f J^/qy/ J - q - Ufy/qJ(0, a^fq) - %\fjy/q\/J - qJ(0, ay/q)), 

nf°) J = 2V7(<Z 3/2 J(l, a^) 2 a 2 ^J - q + Ay/qy/ J - q + Ay/qy/j - qJ(0, a^f 
+4J(1, ay/q)aqy/j - qJ(0, ay/q) - 47(1, ayfq)aqyj J - q 
-&y/qy/J^q~J{Q, ay/q)) + /(a 2 ^ 2 J(0, a^f - a 2 q^ 2 J(0, ay/q) 
+2aq 2 J(l, ay/q) - 2 J(l, a^)aq 2 J(0, a^/g) + 2Ja 2 q 3/2 J(0, ay/q) 
-2Ja 2 q 3/2 J(Q, a^f - 4JaqJ(l,a,y/q) + 47.7(1, aVg)aq J(0, a^)), 

nf) J = (/a 2 g 5 / 2 J(0, a^) 2 - /«V /2 J(0, aV?) - 2fJ(l,ay/q)aq 2 J(0, a,/?) 

+2/ag 2 J(l, aV?) + 2v / Jg 3/2 J(l, ay /q) 2 a 2 y/ J - q + 2 J/a 2 g 3/2 J(0, a^g) 
-2 Jfa 2 g 3/2 J(0, a^) 2 + 4J/J(1, a^)ag J(0, a^) - AJfaqJ(l,ay/q) 
+8Vl^/q^J- qJ(0, a^/q) 2 - 8\/~JJ(l,ay/q)aq\/ J — q + 8^/1^/q y/j - q 
+8VJJ(1, a^/q)aqy/j - qj{0, ay/q) - l&yfjy/qy/j - qj(0, ay/q)), 

n (23)J = fa 2 q5 /2 J(Q; _ /a 2 g 5/2 J(0j ^ + 2/J(1 , a ^ )ag 2 

-2/ag 2 J(l, ay/q) J(0, a^g) + 2\/jg 3/2 J(l, ay/q) 2 a 2 y/T~q 
+2Jfa 2 q 3/2 J{Q,ay/q) - 2Jfa 2 q 3/2 J{0,a^) 2 
+4JfaqJ{l,ay/q)J(0, ay/q) - AJfaqJ(l, ay/q) 
-8\/jJ(l, ay/q)aqyfj - q + 8VJy/qy/j - qj(0, a^g) 2 
+8VjJ(l,ay/q)aqy/j - qJ(Q,a^/q) + %\/j y/qy/ J - q 

-Vb\/lyfq\/J - qJ(0, ay/q). (34) 

From Figs. ([7]) and (|SJ), we see how the pure components of the metric tensor behave, when they are plotted 
against q and the angular momentum. For the choice L = 22500, p = 5/6, 6 = 1 and a = 100000, the gg-heat 
capacity as depicted in the Fig. (JTJ) shows a large variation of the minima and maxima, which respectively take 
in an order of the amplitude +10 10 and — 10 11 and occur for q,J£ (0,4). In this case, we further observe for 
all J £ (0, 4) that the underlying quarkonia become highly unstable as q positively approaches the origin. 

On the other hand, the J J- heat capacity in Fig.© shows only negative variation of order — 10 5 , which 
occurs when q, J tend to the origin. In this case, we notice that all local interactions are present in the region 
q £ (0,4) and J £ (—4,4). The strength of these local interactions depends on the domain chosen in the (q, J) 
space. These fluctuations signify the thermodynamical interactions present in the strongly coupled quarkonia. 

The corresponding plot for the determinant of the metric tensor is depicted in Fig. ([9]). For the above choice 
of the parameters, the Fig.© shows the graphical nature of the determinant of the metric tensor. These plots 
explicate the thermodynamically (un)stability regions for the underlying massive quarkonia lying in a given 
QCD phase. Combining the effects of all fluctuations of the {q, J}, we observe that the quarkonia are stable 
for q := Q 2 £ (1,4). In general, the global stability requires that the determinant of the metric tensor must be 
positive definite, which in the present case transform as 

n J g := n^ J + n^ J l[p) + n^ J l[pf < 0. (35) 

In this case, it turns out that the thermodynamic curvature may be written as the series of the charmo- 
nium logarithmic factor l(p) and Bloch-Nordsieck logarithmic factor f(q, J) of rotation as the coefficient of the 
expansion. Systematically, the exact expression for the scalar curvature takes the form 

^) = ^^E B « X (^))"> (36) 
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where the B n in the numerator of the scalar curvature are polynomials in p, whose coefficients are the functions 
of the Bloch-Nordsieck logarithmic factor f(q,J). While, the denominator of the scalar curvature precisely 
takes the numerator of the determinant of the metric tensor as its square. The quantitative properties of the 
scalar curvature and the Riemann curvature tensor remain similar, as we shall discuss in the next subsection. 

4.2 Quarkonia in QM-plane 

In the present section, we analyze the nature of the fluctuating quarkonia generated by the momentum Q, 
and mass m. To focus on the general case, we choose the variable Q as the transverse momentum with the 
understanding that k = k± and the mass considered as an arbitrary real parameter of the system. Following 
the convention of Bloch-Nordsieck rcsummation, the strong QCD coupling can be expressed as 



a,, n 1 P i ,m + \/m 2 — k 2 ... . . 

b ln(l + p(k 2 /L)P) m - Vm 2 - k 2 

When the Bloch-Nordsieck resummed strong QCD coupling A(k, m) is allowed to fluctuate as a function of 
the {fc, m}, we may again exploit the definition of the Hessian function Hess(A(k, m)). Herewith, we find that 
the components of the metric tensor are given by 

,. - E (n {0)M +n {1)M l(v) + n {2)M l(v) 2 ) 

9kk ~ W(p) 3 exp(2G(p))fc6(m2-fc2)3/2^" + n " 'W +n n «W h 

, - ?? (n {0)M + n {1)M l(v)) 

9km " bl(p) 2 exp(l(p))k(m 2 -k 2 r/ 2[l2 + 12 m) ' 

**» - mB^kW^ - 1+J{0 ^ (38) 

where the coefficients {ji"p M 1 n^ M } appearing in the determinant of the metric tensor factorize as 



(1) M 
l ll 

(2) M 
'11 



2(fe7i)VKi 



'2n[\ 3)M p(k 2 /Lr), 

2 (22)M, 



+ 2pn[ 2 1 1)M (k 2 /Lr + p 2 n\r> lu (k 2 /L) 



\2p 



(39) 
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As it might be expected, we find without any approximation that the factors in the numerator of the local 
pair correlation associated with the fcfc-components can, as before, be easily expressed as 



it 



(° )M = 8p 4 f((k 2 /L) p ) Vm 2 - fc 2 (m 2 - k 2 + J(0, ak)k 2 - J(0, ak)m 2 ), 

™ )M = fVm 2 - fc 2 J(0, ak)k 2 - f^/m 2 - fc 2 J(0, ak)m 2 + Am 3 - 4m 3 J(0, ak) 



+Amk 2 J(0,ak) + f\Jm 2 - k 2 m 2 - j 'y 'm 2 - k 2 k 2 - Amk 2 
-2/ J(l, ak)ak^m 2 - k 2 m 2 + 2/ J(l, ak)ak 3 y/m 2 - k 2 , 
i[\ 3)M = fVm 2 - k 2 k 2 + f{{k 2 /L) p ) 2 Vm 2 - k 2 m 2 - f ((fc 2 / Lff V m 2 - k 2 k 2 
-f^m 2 - k 2 J(0, ak)k 2 - fyjm 2 - k 2 m 2 + fyjm 2 - k 2 J(0, ak)m 2 
-f((k 2 /L) p ) 2 ^m 2 - fc 2 J(0,afc)m 2 + f((k 2 /L) p ) 2 ^m 2 - fc 2 J(0,afc)fc 2 



ii 



-2fJ(l,ak)a((k 2 /L) p ) 2 kVm 2 - k 2 m 2 + 2/ J(l, afc)a((fc 2 /L) p ) 2 fcVm 2 - k 2 
+Amk 2 ((k 2 /L) p ) 2 J(0, ak) - Am 3 {{k 2 / L) p ) 2 J(0, ak) + 4m 3 ((fc 2 / L) p ) 2 
~Amk 2 ((k 2 /L) p ) 2 , 

(20)M = f ak 3^ m 2 Z k 2 J{l,ak) - fak^/m 2 - fc 2 J(l,afc)m 2 + fa 2 k 2 ^/m 2 - fc 2 J(0,afc)m 2 

-fa 2 k 4 \/m 2 - k 2 J(0, ak) + 2m 3 + Amk 2 J(0, ak) - Am 3 k J(l, ak)a - Amk 2 

~2m 3 l(p) 2 J(0, ak) + Amk 3 J(l, ak)a, 
n { £ )M = fak 3 y/m 2 - k 2 J(l,ak) - fa 2 k 4 y/m 2 - k 2 J(0, afc) - fak^/m 2 - fc 2 J(l, ak)m 2 

+fa 2 k 2 ^m 2 - k 2 J(0, ak)m 2 - 2m 3 J(0, afc) + 2to 3 - 4m 3 fc J(l, ak)ap - 4rafc 2 

+4mfc 3 J(l, ak)a + Amk 2 J(0, afc), 
n (22)A/ = j- ak 3^ m 2 _ k 2 J{l,ak) ~ fa 2 k 4 ^/m 2 - fc 2 J(0,afc) + fa 2 k 2 ^/m 2 - fc 2 J(0,afc)m 2 



-fakyjm 2 - k 2 J(l, ak)m 2 + 2m 3 - 2m 3 J(0, afc) - 4m 3 fc J(l, afc)a 

-4mfc 2 + 4mfc 2 J(0, afc) + 4mfc 3 J(l, afc)a. (40) 
In the present case, the factors of the numerator of the fcm-component of the metric tensor are 
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n^ M = 2p 2 (k 2 /L) p (k 2 + J(0,ak)m 2 - J(0,ak)k 2 -m 2 ), 
n { i2 M = k 2 -k 2 J(0,ak) + J{l,ak)akm 2 - J{l,ak)ak 3 
+J(1, ak)akm 2 p(k 2 / L) p - J(l, ak)ak 3 p{k 2 / L) p 

-k 2 J(0, ak)p(k 2 /L) p + k 2 p(k 2 /L) p . (41) 

When the fluctuations of the quarkonia are described in terms of the transverse momentum fc and mass 
m, we observe that the principle components of the metric tensor, signifying self pair correlations, are positive 
definite functions in a non-trivial range of k. For the configurations with the same sign of the index p and 
constant &, it turns out that the local stability requires that (i) qq- fluctuations satisfy the constraint 

n i T + n^ M l(p) + n^ M l( P r>0, (42) 

and (ii) mm- fluctuations pose a constraint on the Bcsscl function such that the underlying transverse momentum 
of the configuration must be limited to the condition of J(0, ak) > 1. 

The determinant of the metric tensor turns out to be in a similar rational form, as it was obtained in 
the previous subsection. However, the simplifications are relatively straightforward, which in turn implies the 
following compact expression for the determinant of the metric tensor 

» - -s^^v^'''?" + «r <« + »r "<&>«. (43) 

with the following factorization 

»r = 2 P 2 (fc 2 /i) p (K i2)M +p(( fc2 /i) p )4 i3)M )' 

nf M = nf» M + 2nf^ M p{k 2 /LY + nf 2 ^ 1 ((k 2 / Lffp 2 . (44) 
In this case, it is relatively easy to obtain that the coefficients are given by 



n 9 



= 8p 4 (v / m 2 - k 2 ({k 2 /L)Pf + ^m? - k 2 {(k 2 /Lf) 2 J(0, akf 
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Figure 12: The determinant of the metric tensor plotted as a function of the transverse momentum and mass 
k, to, describing the fluctuations in massive rotating quarkonia. 



-2mf({k 2 /L)P) 2 J(0, ak) + mf((k 2 /L) p ) 2 + mf((k 2 /L) p ) 2 J{0, akf 
-2^m 2 - k 2 ((k 2 /L) p ) 2 J(0, ak)), 
mf - 2mpf + 2to/J(1, ak)akJ(0, ak) + mfJ(0, ak) 2 
+4mpfJ(0, ak) - 2mpfJ(0, ak) 2 - 2m/ J(0, ak) 



+4\/m 2 - k 2 J(0, ak) 2 + A^Jm 2 - k 2 
+4y / m 2 - k 2 J(0, ak)J(l, ak)ak - 2mfJ(l,ak)ak 
-A^m 2 - k 2 J(l,ak)ak - 8\/m 2 - k 2 J(0, afc), 
4 13)M = to/ + 2to/ J(l, afc)afcj(0, ak) - 8\/m 2 - fc 2 J(0, afc) 



+4-\/m 2 - k 2 + 4\/to 2 - k 2 J(0, afc) J(l, afc)afc 

+4\/to 2 - fc 2 J(0, afc) 2 - 4a/to 2 - fc 2 J(l, afc)afc 
-2m/J(0, afc) - 2mfJ(l,ak)ak + m/J(0, afc) 2 , 
nf 0)M = m(fa 2 k 2 J(0,ak) ~ fak,J(l,ak) + fak,J(l,ak)J(0,ak) 



-/a 2 fc 2 J(0, afc) 2 ) + 2-^/m 2 - fc 2 (l + J(0, afc) 2 
+ J(l, afc) 2 a 2 fc 2 + 2fc J(0, afc) J(l, afc)a - 2fc J(l, afc)a - 2J(0, afc)), 
?i^ 21)M = m(fa 2 k 2 J(0,ak) - f a 2 k 2 J (0 , ak) 2 + /afc J(l, afc) J(0, afc) 



-/afc J(l, afc)) + 2v/to 2 - fc 2 (l + J(l, afc) 2 a 2 fc 2 
+ J(0, afc) 2 - 2fc J(l, afc)a - 2 J(0, afc) + 2fc J(l, afc)aj(0, afc)), 
t^ 22 ) m = TO/(afcJ(l,afc)J(0,afc) +a 2 fc 2 J(0,afc) -afcj(l,afc) 



-a 2 fc 2 J(0, afc) 2 ) + 2v/to 2 - fc 2 (l + J(l, afc) 2 a 2 fc 2 

-2 J(0, afc) + 2fc J(l, afc)a J(0, afc) + 2 J(0, afc) 2 - 2fc J(l, afc)a). (45) 

As per the requirement of the positive dcfiniteness of the determinant of the metric tensor, the global stability 
of underlying quarkonia with fluctuating {fc,TO,} leads to the following constraint 

n^ M + n^ M l(p)+n^ M l(p) 2 < 0. (46) 
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The globally covariant Riemann curvature tensor Rkmkm offers a similar property, as described in the foregoing 
subsection. Specifically, we find that the Rkmkm has various factors of {(m 2 -fc 2 )"/ 2 |n G Z} and various powers 
of the logarithmic coupling l(p). Some of the interesting terms of Rkmkm are -600m 6 (m 2 - /c 2 ) n / 2 ;(p) 4 fc 18 , 
+712m 18 (m 2 - fc 2 ) 9 / 2 Z(p) 4 fc 8 , ... and -1244m 28 (m 2 - k 2 ) 7 / 2 l(p) 4 k 8 . It is worth mentioning that R km km 
contains approximately 50,000 terms whose presentation is quite lengthy for the present paper. Nevertheless, 
the globally invariant scalar curvature and globally covariant Riemann curvature tensor can straightforwardly be 
obtained. As described in the previous subsection, the globally covariant physical properties of two parameter 
quarkonia follow directly from the product of two geometrically invariant quantities, viz. R and |j<?||. As per 
our consideration of the foregoing subsections for the globally invariant scalar curvature R and determinant of 
the metric tensor ||<7||, the covariant properties, viz., Rkmkm follow directly for the two parameter quarkonia. 

From Figs. ([TU|) and (|TT|) . we see that the graphical views of the heat capacities arc defined as pure components 
of the metric tensor. The plots have been depicted against the transverse momentum k and mass m. In the 
present case, the choice of parameters is made as L = 22500, p = 5/6, b — 1 and a = 100000. Notice that the 
graphical nature of all thermodynamic interactions remains the same, as long as the parameter b has a fixed 
sign. 

The fcfc-heat capacity Fig. (|10[) shows a large variation of minima when going from a negative k for a given 
positive mass m £ (0, 5). The order of instability turns out to be as high as — 10 16 . In this case, we further find 
in the other quarqants that the underlying quarkonia remains nearly stable in three distinct pathes. 

On the other hand, the mm-heat capacity Fig. ([lip indicates a similar behavior of the negative variation. In 
contrast to the foregoing case, in this case we find that the local interactions are present only in four disjoint 
regions, all of which exclude the origin. As expected, the strength of the above interactions depends on the 
domain chosen in the (fc, m) plane, and thus the nature of the stability of quarkonia with the fluctuating {k, m}. 

The corresponding plot for the determinant of the metric tensor depicted in the Fig.(fl"2|). shows the graphical 
nature of the determinant of the metric tensor. These plots explicate the of regions of the thermodynamic 
(un)stability for the underlying massive quarkonia lying in a given QCD phase. Specifically, the significance 
of the fluctuations of {fc, m} is that they show thcrmodynamical interactions in the strongly coupled massive 
quarkonia. The above figures indicate that the quarkonia can be highly unstable even in the linear Regge regime. 
As the non-linear effects become stronger and stronger, it is thus expected that the thermodynamic instability 
and correlations would grow further. This motivates us to extend our analysis to the general quarkonium 
configuration. 

4.3 Generic Quarkonia 

In the present subsection, we analyze the properties of general massive rotating quarkonia, when all parameters 
of the theory are allowed to fluctuate. To do so, let us consider the scale q, index of effective potential p 
and angular momentum J as the parameters of the present interest. In the framework of the Bloch-Nordsieck 
resummation, the strong QCD coupling takes the following form 

After some simplification, we find that the components of the metric tensor are 



9qp 



ibl(p) 3 exp (2l{p))q b / 2 (J - qf' 2 

I ( n ^ G + n ^ G l( v ) + n {2)G l(v) 2 

2W(p)3 exp (2l(p)) q z/ 2 (J - q y/2 [ni2 + " 12 l[ P) + n u HP) 



a . - E (n i0)G +n (1)G l(v)) 

9qJ ~ 2W(p) 2 exp(/(p))g3/ 2 (J- g )3/ 2 Ji/2^i3 l[P)) > 

ftp = Mp^^ M + AO.^^ + nS^)), 

^ ^ - blip) 2 exp (l(p))(J - WW + m aV ~ q))(ng)G + mn ™ G) 



gjj 



2bl{p)(J-q)3/ 2 J3/- 



:(-l + J(0,ov/(?)))(2J-g), (48) 
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where the coefficients {n^ G , n^ G ', n^ G , n^ } appearing in the components of the metric tensor factorize 
as mentioned before. In the powers of p, the exact factorizations are given as follows: (i) the gg-componcnt 

= n^ G + 2 P { q /Lfn^ G +P \{ q /Lf)\ (49) 

(ii) the gp-component 

n[T = p{q ,LYn^ G +P \ {q /Lyfn^ G , 

= n^ G + 2 P { q /LYn^ G +P \{ q /LYfn^ G . (50) 

Without any approximation, we obtain that the factors in the numerator of the qq-components are given by 

<i )G = 8p 4 f(( q /L)n 2 ^J~ q ( q 3/2 + VvJ(0,aVv)J-VvJ 
- q 3 ^J(0,a^)), 

n[\ 2)G = }y/J~ qq 3/2 - }^J~ qq 3/2 J{0,a^)- fy/J^^J 
+fJ{l 1 a^/ q )a q ^f7~ q J - fJ(l,ay/q)aq 2 y/T~q~ 

+fVj Zr q~Vv J (Q, a^)J - 2VJ q 3/2 J(0, a^) 
-2 J 3 / 2 ^ + 2 J^ 2 ^J(Q, ay/q) + 2%/jV 3 / 2 

+pfV J - <iVq J - pfV J - iV^ J (°i a V?) J 

-p.f^J~qq 3/2 + P fVj~ qq 3/2 J{0,a^f q ) 1 
= /J(l, a,yfq)aq\J J - q J - fJ(l, a,y/q)aq 2 \J J- q 

-fx/T^y/qJ - f^J~ qq 3/2 J{Q, aV?) 

+.f^J^^J(0, a^)J + f^J~ qq 3/2 

-2 V7g 3/2 J(0, ay/q) + 2J 3/2 ^J(0, ay/q) 

+2VJ q 3/2 - 2J 3 / 2 V9, 
rif? )G = fa(aq 5 / 2 y/j - qj(0, a^) - aq 3 ' 2 y/j - qj(0, a^)J 

+2 qs /j-qJ(l,ay/q)J - 2 q 2 ^ J - qj{l, ay/q)) 

+2V7(2 Jy/qJ{0, ay/q) - 2J(l,ay/q)aq 2 

+2 JJ(1, ay/q)aq + 3q 3/2 - 4 J 3/2 ^/q~ - 3q 3/2 J{0, ay/q)), 
nii 1)G = fa 2 q 5 ' 2 JJ^q~J{Q,ay/q)- fa 2 qV 2 ^fj^q~J{U,ay/q)J 

+2faqy/j - qj(l, a^)J - 2faq 2 y/ J — qJ(l, ay/q) 

-4J 3/2 y/q~ + 4J 3/2 J(l, a^/q)aq - iVJj(l, a^acj 2 

+4J 3 / 2 ^J(0, ay/q) + 6VJ q 3/2 - 6VJ q 3/2 J(0, ay/q), 
n (22)G = fa 2 q 5 / 2 y/j^J(0,a^)-fa 2 q 3 ^ 2 ^J^ 1 J(0,a^)J 

+2/J(l, ay/fiaqy/j -qJ- 2fJ(l,a^ q )a q 2 y/j - q 

-4VJj{l,a^)aq 2 + 4J 3/2 ^J(0, a^) 

-AJ 3/2 y/q + 4J 3/2 J(l, a^)aq + 6\fj q 3/2 

-6\/jg 3/2 J(0,aVg). (51) 
In a parallel way, the factors in the numerator of the gp-components turn out to be 

< 2 )G = 4p 3 f(( q /L)n 2 (x/j~qVv-Vj^VvJ(^*Vv) 
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i\fT~q^Jqhi(q/L) - p\J J - qyfqj (0, ay/q) ln(q/L)), 



,(H)G 



(12)G 



= 2^7 y/q - 2p 2 fy/j^q~y/q]n(q/L) - &p.f^J~q^q 
+2p 2 f^J~qy/T l J{Q, a% /q) Hq/L) - 2VJy/qJ(0, a% /q) 
-fJ(l, ay/q)aqx/J - q + 2pVJy/qln(q/L) 
-2pVJy/qJ{0, ay/q) ln(g/L) + 6pfy/J=qy/qJ{0, a^q) 
—pfJ (1) a x /q)aqy / J — q ln(q/ L), 

= 2y/jy/q- / J(l, ay/q)aqy/J -q + 2pVJ ^ \n(q / L) 
-2V1^J(0, ay/q) - 2pV7^J(0, a^) In(g/L) 

+*pfy/J -qVqJ(0,a\fi) - Ipfy/J-lVv 
-pfJ (1, a^/q)aqy/j - q]n(q/L), 

= faqy/J - qJ(l,ay/q) + 2yfjy/qJ{Q, ay/q) - 2\/jy/q, 



''12 

n i2°^ G — ag-\/ ^ — q J (I, ay/q) + 2\fly/qJ{0 1 ay/q) - 2^/ly/q, 

(22) G 
"12 



2yfjy/qJ%ay/q)-2y/jy/q + faqs/J^lj(l,ay/q). (52) 
The factors in the numerator of the g ./-components are expressed as 

<3 )G = 2p 2 (g/L)f (g 3 / 2 + ^J(0, a^g) J - g 3 / 2 J(0, a^/q) - y/qJ), 

"is G = 9(9 1/2 -g 1/2 J(0,aV?) + J(l,aV«)oJ-J(l,aVg)a«) 

+p( g /L) p (<Z 3 / 2 - g 3/2 J(0, a^) + J(l,ay/q)aqJ - J(l, a^)ag 2 ). (53) 

The corresponding factors in the numerator of the pp-componcnts are given by 

n 2 ° )G = -2p( g /L) 2 P(l + 2pln(g/L)+p 2 ln(< 7 /i) 2 ), 
4f - (g/L)P(2 + 4pln( (7 /L)+p 2 ln(g/L) 2 ) 

+p( q /LfP(l + 2pln(q/L)). (54) 

Finally, the factors in the numerator of the p ./-components are 

4T = - P (q/Lr-p 2 \n(q/L)(q/Lr, 

ng )G = l+p(q/Ly. (55) 

Herewith, wc observe that the principle components of the metric tensor, signifying self pair correlations, 
remain positive definite functions in non-trivial intervals of the parameters. The local stability of the configu- 
ration requires the following three simultaneously constraints. For the same sign of {p, b, the thermodynamic 
stability enforces: (i) the qq- fluctuations satisfy 

n^ G + n^ G l(p)+n^ G l(p) 2 , (56) 
(ii) the J J- fluctuations remain within the limiting values of the Bcssel function 

J(0,ay/q) > 1, 2 J < q, 

< 1, 2J>q, (57) 

and (iii) the pp- fluctuations satisfy 

4f +4fi(?) > 0, J{0,ay/q)>l, 

< 0, J(0,ay/q)<l (58) 

for the same sign of {/,&}. Correspondingly, the three dimensional views of the heat capacities are shown 
in the following plots for the general rotating massive case. For the purpose of the present and subsequent 
diagrammatic views, wc shall focus on the case of p = 5/6 and do not take any approximation against the 
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Figure 13: The gg-component of the metric tensor plotted as a function of q and the angular momentum, J, 
describing the nature of the heat capacity of the JJ-component for general massive rotating quarkonia. 



Bessel function. Thereby, we graphically illustrate the local stability of the rotating quarkonia, which is in 
accordance with the above positivity constraints of the specific heat capacities. 

As a function of scale parameter q := Q 2 and angular momentum J, the gg-component is shown in the 
Fig.([T3"l). As mentioned before, let us illustrate the present case for the same choice of the parameters, viz., L = 
22500, 6 = 1 and a = 100000. Consequently, we find that the graphical nature of the underlying thermodynamic 
nature of the quarkonia remains the same, as long as the parameter b has a fixed sign. 

The pp-hcat capacity is depicted in the Fig.(fT4|). For a small g, this shows a large minimum of depth 10 6 , 
in the regime of the angular momentum J G (0,5). It turns out that the underlying pp-fluctuations of the 
quarkonia remain nearly stable in the other quadrants. 

Similarly, the JJ-heat capacity is shown in the corresponding Fig. (fT5")) , which indicates a similar behavior 
of the negative variation, as mentioned in the case of the Q J-plan. The silent feature of the JJ-fluctuations is 
that we find the two distinct local behaviors for J > and J < 0. In the limit of vanishing q and J, the local 
interactions reach an order of 10 4 . For fluctuating {q,p, J}, the strength of the thermodynamic interactions 
depends on the domain chosen in the (g, p, J) plane, and thus so does the thermodynamic stability of the 
quarkonia. 

We further observe that the Gaussian fluctuations of the most general quarkonia comply with expected 
thermodynamic stability. Physically, the heat capacities under the fluctuations, which are defined as the self- 
pair correlations, remain positive quantities in the domain of the parameters and thus form the well-defined 
basis of the manifold (M3, g). Subsequently, our computation demonstrates the intrinsic geometric nature of the 
parametric fluctuations. The quarkonia remain globally stable under the fluctuations, if the associated principle 
minors, viz., {p2,Ps\ remain positive definite functions. Herewith, the gp-surface is stable if we have a positive 
surface minor 

»" - -im^wm¥w^ (n "' a + "» + '° i(f,)2 + "™ G,< '' )3, • (59) 

where the rig has the following illuminating factorization 

4 0)G = 8/ 2 (( 9 /l)w4 01)g +M 02)G + P 6 4 03)G )- (60) 

(01)G (02)G (03)G-| 



Interestingly, the factors {n s >n-s , tIq } can be expressed as 
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Figure 14: The pp-component of the metric tensor plotted as a function of the scale q, and the angular momen- 
tum, J, describing the nature of the heat capacity of the JJ-component for general massive rotating quarkonia. 




q 3 / 2 ^J~q + q 3 ' 2 ^T^J(0, a^) 2 - J^^J^l 
-J^^I~qJ(Q, a^) 2 - 2q 3 l 2 ^J~qJ{Q, a^q) 
+2Jy/qy/J-qJ(0,ay/q), 

2q 3/2 ^/J -qln(q/L) - 2J^J-qJ(0, a^) 2 \n(q/L) 
+2q 3 / 2 ^J -qJ(0, a^) 2 ln{q/L) + 3q 3/2 y / J-q 
-3y/qy/ J - qJJ(0, a^f + Zq 3 ' 2 y/1 - qj(0, a^) 2 
+4J^y/7^J(0, aV?) Hl/L) - 2JJqy/j - qln(q/L) 
-6q 3 / 2 ^/j-qJ(0, a^q) + 6y/qy/j-qJJ(0, a^q) 
-2,^^J~qJ - 4g 3 / 2 ^J~qJ(Q, ajq) ln(g/L)), 
q 3/2 ^J~qHq/L) 2 + q 3 ' 2 ^T~qJ{0, a^/q) 2 Hq/L) 2 
-J^^J~qJ{Q, a^) 2 Hq/L) 2 - 2^^J~qJ Hq/L) 
+2q 3 ' 2 ^J~q,J(Q, a^/q) 2 \n(q/L) + 2J^^I~qJ{Q, a^) Hq/L) 2 
+2q 3 ' 2 ^J -qHt/L) - 2q 3 / 2 ^J-qJ(0,a^) Hq/L) 2 
-2^^T~qJJ{0, a^) 2 Hq/L) + ^^J~qJJ(Q, ay/q) ln(g/L) 
-J^^I~qHq/Lf -4g 3 / 2 v/7^J(0, a^)Hq/L). (61) 

It is worth mentioning that the factors of the /(p)-terms can further be expressed as 



« G = {q2/LfVnf )G + {{k 2 /L) 3 ?nf )G ) (62) 
with the following additional structures 



n (X2)G = n (120)G + n (121 )Gf + n ^2) G/2j 

nf) G = 4 13 °) G + 4 131 ) G / + 4 132 ) G / 2 . (63) 
We find in particular that the 120-factor reduces to the following explicit expression 



„(01)G _ 
n s — 



4 G = MVq^J~qJp 2 + VqVJ^qp i JJ(o,a^) 2 in( q /L) 2 

-2q 2 ^Jp 3 fJ(l,a^)aHq/L) - 2^J-qJJ(0,a^)p 2 
+2^^J~qp 3 JJ{Q, a^q) 2 Hq/L) - 2^^J~^p 4 JJ(0, a^) ln(g/L) 2 
-^^J~qp 3 JJ{0,a^q)Hq/L)). (64) 

The 121- and 122-factors factorizc as follows 



4 21 > G = p 2 nf 12)G + 2p 3 nf 13)G +/4 1214)G , 

4 122)G = p 2 nf 22)G + 2p 3 4 2123)G + /4 2124)G + 4p 5 4 2125)G (65) 
with the following expressions for the sub-factors 



rig 212 ^ G = Vl(JqaJ(l,a^/q) — q 2 J(l,a^q)a + q 2 J{l,a^/q)aJ{Q,a^/q) 
-JqaJ(l,a^q)J(0, a^) - 2J^ - 2J^J(0, a^) 2 + Zq 3 ' 2 
+'iq 3 ' 2 J(0, ajq) 2 + 4 J^J(0, a^q) - &q 3 ' 2 J(0, a^q)), 

4 2U3)G = Vj(q 2 J(0 1 a^)Hq/L)J(l,a^)a + JqJ{l,a^)aHq/L) 



26 



- 7g7(l, ay/q)aJ(0, ay/q) ln(g/L) - 2J^\n(q/L) + 3g 3/2 
-2 7^7(0, a^g") 2 ln(g/i) - 37^/g" + 3g 3/2 7(0, a^) 2 ln(g/7) 
+3g 3/2 7(0, a^) 2 - 3 7^7(0, a^g") 2 + 3g 3/2 hi(g/Z) 
+4 7^7(0, ay/q) ln(g/Z) + 6 7^7(0, a^g) 
-6g 3 / 2 7(0, a^ff) ln(g/7) - 6g 3 / 2 J(0, a^g")), 
n aai4)G = Jj( q 2j(i^ a ^q)aJ(0, a^g) \n(q/L) 2 + Jq]n(q/L) 2 J(l,a^/q)a 
-q 2 ln(g/L) 2 7(l, ay /q)a - 7g7(0, ay/q) Hq/L) 2 7(1, ay/q)a 
-2 JVgJ(0, a^/g) 2 ln(g/L) 2 + 3q 3/2 7(0, aVg) 2 m(g/L) 2 
+3g 3/2 ln(g/7) 2 + 47^7(0, ay/q) ln(g/L) 2 
+4g 3 / 2 J(0, a^) 2 ln(g/7) + 4g 3 / 2 ln(g/Z) - 4J^ln( g /L) 
-47^7(0, a^) 2 ln(g/L) - 6q 3 ' 2 7(0, a^) ln(g/L) 2 

+8 J^gJ(0, a^g) - 8g 3/2 7(0, a^g) m(g/7)). ( 66 ) 

The 122-factors take five powers of p, whose precise expressions take the form 



n (2l22)G = agv /j _ g(gV2 j a j(l i a ^)2 _ q W a J(l,ay/q) 2 + 27g 1/2 a7(0, ay/q) 2 
-2Jq ( l/2)aJ(0, ay/q) + 2q 3/2 aJ(0, ay/q) - 2g 3/2 aJ(0, a^g) 2 
-4JJ(1, a^g) J(0, a^g") + 4g7(l, a^ff) J(0, a^) - 4g7(l, a^g) 
+477(1, aVg)), 

4 2123) G = VJ^~q(q 3/2 Ja 2 J(l 1 a^q-) 2 ln(q/L) - q 5 / 2 J(l,ay/q) 2 a 2 ln(q/ L) 

+27g 3/ V J(0, ay/q) 2 hx{q/L) + 2g 5/2 a 2 7(0, ay/q) Hq/L) 

-2g 5/2 a 2 7(0, ay/q~) 2 ln(g/L) - 27g 3/2 a 2 7(0, ay/q) Hq/L) 

+4g 2 a7(l, Oy/q)J(0, a^g) ln(g/7) - 47ga7(l, a,y/q)J(0, ay/q) Hq/L) 

+47ga7(l, ay/q) ln(g/L) - 4g 3/2 + 4Vg7p 5 7(0, ay/q) 2 Hq/L) 

— 4g 2 a7(l, ay/q) \n(q/L) + 6g7a7(l, ay/q) — 6g7a7(l, ay/q)J(0, ay/q) 

+6g 2 a7(l, ay/q) 7(0, a^g) - 6g 2 a7(l, ayfi) + 8g 3/2 7(0, a^g) 

-87Vg7(0,aVg)), 

4 2124) G = v^7^(g 3 / 2 77(l,aV5) 2 a 2 ln(g/L) 2 -g 5/2 7(l,aVg) 2 a 2 l n (g/i) 2 

+2g 5 / 2 a 2 7(0, a^/g) ln(g/L) 2 - 2g 5 / 2 a 2 7(0, a^g) 2 ln(g/£) 2 
-27g 3/2 a 2 7(0, a^g) ln(g/L) 2 + 27g 3/2 a 2 7(0, a^) 2 ln(g/L) 2 
-47ga7(l, 0^5)7(0, a^/g) ln(g/L) 2 + UqaJ(l,ay/q) ln(g/L) 2 
+4g 2 a7(l, ay/q) 7(0, a^g) ln(g/£) 2 - 4g 2 a7(l, a^/g) ln(g/L) 2 
+8g 2 7(0, oV?) In(g/L)a7(l, ay/q) - 8g77(0, a^g) In(g/L)a7(l, a^/g) 
+8g7ln(g/L)a7(l,o 1 /g) - 8g 2 7(l, ay/q)aHl/L) 
+16JVg7(0, a% /g) 2 ln(g/L) + 167Vghi(g/7) 
-16g 3/2 7(0, ay/q) 2 ln(g/7) - 16g 3/2 ln(g/L) 
+28^77(0, cVg") 2 - 28g 3 / 2 7(0,aV5) 2 
+28y/qJ - 28g 3/2 - 327Vg7(0, ay/q) ln(g/7) 
+32g 3 / 2 7(0, a^g) ln(g/7) - 56VgJ7(0, a^g) + 56g 3/2 7(0, ay/q)), 
nf 2b)G = y/j^iJy/q-Ha/L) 2 + Jy/qJiO, a^) 2 ln(g/L) 2 - q 3 ' 2 ln(g/L) 2 

-2g 3 / 2 7(0, a^g) 2 ln(g/7) - 27VgJ(0, ay/q) ln(g/i) 2 + 2Vg71n(g/L) 
+2g 3 / 2 7(0, ayft) ln(g/L) 2 - 4^77(0, ay/q) In(g/7) 

+4g 3 / 2 7(0,aV?)ln(g/i))- (67) 
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Further, we take advantage of the following factorizations for the 131 and 132 expressions 



4 131)G = v4 i3i3)G + 8 P 4 4 13 14)G + 4p 5 4 1315)G , 

n^ 2)G = p 3 nf 23)G + 2p i nf 2i)G + p 5 n { s 1325)G . 
Subsequently, the factors of 130 and 131 appear in the powers of the p, and are given by 



±yfq~^/T~qp 3 J(l + p 2 J(0, a^) 2 \n(q/L) 2 + J(0, a^) 2 
+p 2 \n(q/L) 2 + 2p\n(q/L) + 2pJ(0, a^q~f \n(q/L) - 2J(0, ay/q) 
-2p 2 J(0, ay/q) Hq/L) 2 - 4pJ(0, ay/q) \n(q/L)), 
J 3/2 qaJ(l,a^/q~) - q 2 \fjj(l,ay/q)a + q 2 aJ{l, dy/q)\fjJ% ay/q) 
-J 3/2 qJ{l,ay/q)aJ{0, ay/q) - 2J 3/2 ^- 2 J 3/2 ^J(0, ay/q) 2 
+3q 3/2 VJ + 3q 3/2 VJj(0, ay/q) 2 + 4J 3/2 ^J(0, ay/q] 
-6q 3 / 2 VJj(Q,ay/q), 

J 3/2 qJ(l, ay /q)a ln(q/L) + q 2 VJj(l, ay/q~)aJ(0, ay/q) ln(g/i) 
- J 3/2 qJ(l, ay/q~)aJ(0, ay/q) Hq/L) - q 2 v f JJ(l,a y /q~)aln(q/L) 
-2J 3 / 2 y^J{Q,ay^) 2 \n{q/L) - 2J 3/2 y^\n(q/L) + 3q 3/2 VJ Hq/L) 
+3q 3/2 VJ + 3q 3/2 VJj(0, a^) 2 - 3 J 3/2 y/q - 3 J 3/2 ^J(0, ay/q) 2 
+3q 3/2 VJj{0, a^/q-f \n(q/L) + U 3/2 ^J{0, a^) \n(q/L) 
+6 J 3/2 ^J(Q, ay/q) - 6g 3/2 VjJ(0, ay/q) - 6g 3/2 VJJ(0, ay/q) ln(q/L), 
J 3/2 q \n{q/L) 2 J(l, ay/q)a - q 2 \ r JJ{l, a^)a Hq/L) 2 

+q 2 VJj(0, a^) ]n(q/L) 2 J(l, ay/q)a - J 3/ V(0, a^) \n(q/L) 2 J(l,a^)a 
-2 J 3/2 V? H<l/L) 2 - 2J 3 ' 2 ^J{Q, ay/qf Hq/L) 2 + 3<z 3/2 V7ln(g/L) 2 
+3g 3/2 \/j J(0, a^f HQ/L) 2 - 4J 3/2 y/qJ(0, a^f \n(q/L) 

~u 3/2 ^Hq/L) + 4g 3/2 V7 J(0, ay ^) 2 Hq/L) + Aq 3/2 ^Hq/L) 
+4J 3/2 V ^J(0, ay/q) HQ/L) 2 - 6q 3/2 \/jJ(0, ay/q) ln(q/L) 2 
-8q 3 / 2 \/j J(0, ay/q) Hq/L) + 8 J 3/2 ^J(0, ay/q) ln(g/L). 

The factors of 132-powers, appearing as the powers of the p, are 

4 1323)G = q 3 l 2 yfJ~qJJ{\,a^) 2 a 2 -cf>l 2 yfT~qa 2 J{l,a^) 2 

-2 Jq 3/2 a 2 y/j - qJ(0, a^) - 2q ( b/2)a 2 y/j - qJ(0, a^f 
+2q 5 / 2 a 2 y/j-qJ{0,ay/q~) - UqJ(l, ay/q)a\J J - qj(0, ay/q) 
+AJqay/j -qj(l, a^) + 2 Jq 3/2 a 2 yj 3 - qj(0, ay/q) 2 
+4q 2 a y / J - qj(l, ay/q) J(0, ay/q) - Aq 2 ay/ J - qj(l,ay/q), 

+2 J^xJT~~q - 2q 3 l 2 yJ,J-q + 2J yfqyj J - qJ(0, a^f 

-2q 5/2 a 2 yjl~ q J(0, a^f Hq/L) + 2q ( b/2)a 2 y/ J - qJ{Q, ay/q) ln(g/L) 

~2Jq 3/2 a 2 \J J - qj(0, Oy/q) Hq/L) + 2Jq 3/2 a 2 y/j - qj(0, ay/q) 2 Hq/L) 

~2q 3 / 2 y/j-qJ(0 1 a^) 2 - 4g 2 J(l, a^ayj J - qHq/L) 

+4JqJ(l,ay/q)ay/ J - q\n{q/L) - 4 Jy/q yj J - qJ(0, ay/q) 

+Aq 3 ' 2 y/J -qJ(0, a^q) - AJqay/J - qJ(l, ay/q)J(0, a^q) Hq/L) 
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„(130)G 



(1313)G 



(1314)G 

■S " 



(1315)G 



+V J(l, a^fq)asfj - qj(0, ay/q) \n(q/L) - 6q 2 \fJ - qJ(l, ay/Fj)a 
+6q 2 y/ J - qJ{0, a.y/q)J(l, ay/q)a + 6qy/ J — qJJ(l,a^/q)a 
—6qy/ J — qJJ(0, ay/q) J (1, a^/q)a, 
4 1325)G = q*/ 2 ^T~qJJ{l,a^) 2 a 2 \n{q/L) 2 -q 5/2 ^T~qJ{l,a^fa 2 \n{q/L) 2 

+2q 5/2 a 2 y/j- qj(0, a^Jq) ln(q/L) 2 + 2Jq 3/2 a 2 y/J-qJ(0, a^/q) 2 ln(q/L) 2 
-2q 5/2 a 2 y/,J - q.J{0, a^) 2 ln(q/L) 2 - 2 Jq 3/2 a 2 y/j - qJ{Q, ay/q) hx{q/L) 2 
+Aq 2 J(l, a^q)a yj 3 - qj(0, ay/q) ln(q/L) 2 + AJqay/ J - qj(l, a^) ln(q/L) 2 
-AJqJ{l,dy/q)ay/j - qJ(0, ay/q) \n{q/L) 2 - 4g 2 J(l, a x /q)ay / J - q\n(q/L) 2 
-8q 3/2 y/J - qHq/L) - 8q 3/2 y/J - q,J(0,a^) 2 \n{q/L) + SJy/qy/1 - q\n(q/L) 
+8Jy/qy/J - qj{0, a^) 2 In(g/L) + 8qy/J-qJJ(l, a^)a ln(q/L) 
-8qy/J-qJJ(0, a^/q) \n{q/L)J(l, ay/q)a - 8q 2 y/ J — qJ(l,a x /q)a ln(q/L) 
+8q 2 y/J -qJ{0, a^/q) ln(q/L)J(l, a^q)a + \6y/qy/J - qJ 
+ 16y/qy/j-qJJ(0, a^) 2 - 16J^y/j - qJ(0, a^) \n{q/L) 
-16q 3/2 y/j-q + 16q 3/2 x/J - qJ(0, a^) \n(q/L) - I6q 3/2 y/J~ gj(0, ay/qf 
-32y/qy/j - qJJ(0, ay/q) + 32q 3 / 2 yj 3 - qj(0, ay/q). (70) 

It turns out that the factor of the l(p) 2 terms can be expressed as 



nf G = {q/LYnf )G + {q/L) 2 Vnf )G + {{q/L) 3 Vnf )G ) (71) 
with the following factorizations 



(21)G 
S 


= nf 0)G - 


h4 11)G /- 


-nf 2)G f 2 


(22)G 
S 


= nf 0)G - 




-nf 2)G f 


(23)G 
S 


= nf 0)G - 


Vnf 1)G f~ 


-nf 2)G f 



and sub-factorizations of the following specific forms 



4 211)G = 4p4 2111)G + 8 P 2 4 2112)G + 2p 3 4 2113)G , 

nf 2)G = 2p4 2121)G + 2p 2 4 2122)G +p 3 4 2123)G . (73) 
The precise sub-factors of the 21-factorizations are 



4 210)G = -8^yfJ~qJp{\ + J(0, ay/q) 2 + p \n(q/L) - 2.7(0, a^) 
+ P J(0, a^) 2 \n(q/L) - 2pJ(0, ay/q) \n(q/L)), 

n {2111)G = 2J 3 ' 2 ^J{Q,a^q) 2 +2J 3 ' 2 ^-2q 3 l 2 y/jJ{Q,a^q) 2 
-3q 3/2 VJ - 4J 3/2 V9 J(0, a^j) + 6q 3/2 VJj(0, ay/q), 

4 2U2)G = J 3/2 qJ{l, ayfq)a J(0, a^) \n(q/L) + q 2 \/jJ(l,ay/q)a]n(q/L) 
-J 3/2 qJ(l,ay/q)a]n(q/L) - q 2 v r JJ(l,a y /q)aJ{0, ay/q) ln(q/L) 
+2J 3 ' 2 ^J{0, ay/q) 2 ln(g/L) + 2J 3 / 2 y/q\n(q/ L) 
-3<? 3/2 VjJ(0, ay/q) 2 ln(g/L) - 3q 3/2 y/jJ(0, ay/q) 2 
~3q 3 / 2 Vj-3q 3 / 2 Vj\n(q/L) + ZJ 3 ' 2 ^ 
+3 J 3 ' 2 V?J(0, a^) 2 - 4 J 3 / 2 V? J(0, aV?) M°/L) 
-6J 3 ^y/qJ(0,ay/q) + 6q 3/2 \/jJ(0, a^) 
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+6g 3/2 VjJ(0, ay/q) Hq/L), 
n (2U3)G = 2J 3 ^^F 1 ln(q/L) 2 + 2q 2 ^JJ(l ia ^)aln(q/L) 2 

-2J 3 / 2 qJ(l,a^)a\n(q/L) 2 - 2q 2 VjJ(l,ay/q~)aJ(0 1 a^) hi(q/L) 2 
+2J 3/2 qJ(l, ay /q~)aJ(0, a^) ln(q/L) 2 - 3q 3/2 \fj\n{q/L) 2 
~3<j 3/2 \/jJ(0, a^f Hq/L) 2 + 4J 3 / 2 ^\n(q/L) 
-U 3 / 2 ^J(0,a^)\n(q/L) 2 - 4q 3 / 2 VJln(q/L) 
-4q 3 /' 2 VjJ(0, a^) 2 \n{q/L) + U 3/2 ^J(0, ay/q) 2 Hq/L) 
+6q 3/2 VJj{Q, ay/q) Hl/L) 2 + 6g 3/2 V7 J(0, ay/q) ln(g/L) 

+2 J 3/2 ^J(0, a^f HQ/L) 2 - 8J 3 / 2 y/qJ(0, a^) Hq/L). (74) 
For the 212-components, we have the following factors 



n [2121)G = y/J~q'(q 5/2 a 2 J(l,a^) 2 + Jq 3/2 a 2 J(0,a^) + q ( -5/2)a 2 J(0,a^) 2 
-q { 5/2)a 2 J(0, ay/q) - Jq 3/2 a 2 J(0, ay/q~) 2 - q 3/2 J a 2 J(l, ay/q) 2 
+2JqaJ(l, ay/q)J(0, ay/q) - 2q 2 aJ(l,ay/q~)J(0, ay/q) + 2q 2 aJ{l, ay/q) 
-2JqaJ(l,ay/q)), 

n (2122)G = y/J^7 1 (q^ 2 a 2 J(l,a^) 2 H<l/L)-q 3/2 Ja 2 J(l,ay/^) 2 ln(q/L) 
-2q 5 / 2 a 2 J(Q, ay/q~) \n(q/L) - 2Jq 3/2 a 2 J(0, ay/q~f \n(q/L) 
+2q 5/2 a 2 J{0, ay/Fj) 2 Hq/L) - Aq 2 aJ{\, ay/q)J(0, ay/q) Hq/L) 
+4JqaJ(l, ay/q)J(0, ay/q) Hl/ L ) + iqJaJ(l, ay/q)J(0, ay/q) 
+2Jq 3/2 a 2 J{0, ay/q) Hq/L) - UqaJ(l, ay/q) ln(g/Z) 
+4g 2 aJ(l, ay/q) ln(g/L) - QqJaJ(l,ay/q) 
-6q 2 aJ(l, ay/q) J(0, ay/q) + 6q 2 aJ{l,ay/q)), 

n {2123)G = y/I^ 1 (q^ 2 a 2 J(0,a^) 2 H<l/L) 2 -q 5/2 a 2 J(0,a^)Hl/L) 2 
+Jq 3/2 a 2 J(0, ay/q) Hq/ l ) 2 - Jq 3/2 a 2 J(0, ay/q) 2 Hq/L) 2 
~2q 2 aJ(l, ayff)J{0, ay/q) ln(q/L) 2 + 2 JqaJ(l, ay/q) J(0, a^) ln(q/L) 2 
+2q 2 aJ(l, ay/q~) Hq/L) 2 - 2JqaJ(l, ay/q) ln(q/L) 2 
+4<? 2 aJ(l, ay/q~) \n(q/L) + iqJaJ(l,ay/q)J(0, ay/q) Hl/ L ) 

-AqJaJ(l, ay/q) Hl/L) - 4g 2 aJ(l, ay/q) J(0, ay/q) Hl/ L )). (75) 
Similarly, it follows that the 22- factors obey the following particular relations 



(221) G „ 2 (2121)G . , 3 (2122)G . „ 4 (2213)G 

n s — 2p n s + 4p n s + 2p n s , 

(222) G 2 (2221)G , „ 3 (2222)G . 4 (2223)G lna ^ 

n s = p n s +2p n s +p n s . (76) 

In the above relations, the factors of the 220 and 221-components are given by 



n (220)G = l6J ^j _ 9 /(2J(0, ay/q) + 2pJ(0, ay/q) ln(q/L) 
-pJ(Q,ay/q-) 2 ln(q/L) -p\n{q/L) - J(0,a^) 2 - 1), 

n [ g 121)G = \/j{2JqJ(l 1 ay/q~)a-2q 2 J(l,ay/q~)a + 2q 2 J(l 1 ay/q~)aJ(0,ay/q) 
-2JqJ(l, ay/q-)aJ(0, ay/q) + 6 Jy/qJ(0, ay/q) 2 - 9g 3/2 J(0, ay/q) 2 
+6Jy/q - §q 3 ' 2 - UJy/qJ(0, ay/q) + 18g 3 / 2 J(0, ay/q)), 

n [2122)G = \/j(2q 2 J(l,ay/q)aHq/L)-2JqJ(l 1 ay/^)aHq/L)-40Jy/q~J(0,ay/q) 
+2JqJ(0, ay/q) ln(g/£)J(l, ay/q)a - 2q 2 J(0, ay/q) ln(g/L) J(l, ay/q)a 
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+12 J^j\n(q/L) + 12Jy/qJ(0, a^f ln(q/L) - 18q 3/2 J(0, a^g) 2 hi(g/L) 
-18(7 3/2 ln(g/i) + 20JVg - 20q 3/2 + 20Jy/qJ(0, a^f - 20q 3/2 J(0, a^) 2 
-24J^J(0, oV5) ln(?/i) + 36g 3/2 J(0, aVs) H a / L ) + 40g 3/2 J(0, a^/?)), 
4 2213)G = VJ{2q 2 ln(g/L) 2 J(l, a^)a + 2J^ln{q/L) 2 - 4<7 3/2 ln(q/L) 
-2Jqln(q/L) 2 J{l, dy/q)d + 2JqJ{0, dy/q) \n{q/L) 2 J(l, a^)a 
-2g 2 J(0, aVs) ln(g/L) 2 J(l, a^a + 2 J^J(0, a^) 2 ln(g/L) 2 
-3g 3/2 ln(g/L) 2 - 3<? 3/2 J(0, a^) 2 \n(q/L) 2 + My/q]n(q/L) 
-4 J^g J(0, aV?) ln(g/i) 2 + 4 J^J(0, a^) 2 ]n(q/L) 
-4q 3 / 2 J(0, a^) 2 ln(g/L) + 6g 3 / 2 J(0, a^) m(<z/£) 2 

-8 J(0, aV?) ln (<?/i) + 8q 3/2 J(0, dy/q) (77) 
In this case, we further obtain that the 222-factors are 



n (222i)G = aq(3q 3 ^ 2 a\JJ - qJ(0, dy/q) 2 + 3Jq 1/2 dy/J - qJ(0, dy/q) 
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-3q 3/2 a\/j- qJ{Q,ay/q) - 3Jq 1/2 a\J J - qj(0,ay/q) 
+Aq 3 > 2 y/J - qj(l, ay/q) 2 a - 4q 1/2 ^ J - qJJ(l, dy/q) 2 d 
-6q\J J - q J {1, dy/q) J {Q, dy/q) - 6J\/ J - qj(l,a.y/q) 

+6Jy/J-qJ{l, dy/q)J(0, dy/q) + Qqy/j-qJ{\, dy/q)), 
n [2222)G = ay/J^q~{2q 5 / 2 J{l 1 dy/q) 2 a ln(g/L) - 2q 3 / 2 JJ(l 1 (Xy/qfa Hq/L) 
+3g ( 5/2)aJ(0, ay/qj 2 \n(q/L) - 3q 5/2 aJ(0, dy/q) ln(g/L) 
-3 Jq 3/2 aJ(0, a^) 2 ]n(q/L) + 3Jq 3/2 aJ(0 1 dy/q) \n(q/L) 
+6JqJ(l, ay/q)J(0, a^) Hq/L) - 6q 2 J(l, dy/q) J(0, aV?) ln(?/L) 
-6JgJ(l, dy/q) ln(q/L) + 6q 2 J{l, dy/q) Hq/L) - WqJJ{l, dy/q) 
+10g 2 J(l, a^) + 10gJJ(l, a^) J(0, aV?) - 10g 2 J(l, a^)J(0, a^)), 
4 2223)G = ay/T=qtf' 2 aJ(0, ay/^f Hq/L) 2 - q 5/2 aj(0, a^) ln(g/L) 2 
+ Jq 3/2 a J(0, a^/s) ln(<?/L) 2 - Jq 3/2 aJ(0 1 dy/q) 2 Hq/L) 2 
+2q 2 J(l, aV«) ln(g/i) 2 - 2g 2 J(l, a^) J(0, a^g) ln^/L) 2 
-2JqJ(l, dy/q) Hq/L) 2 + 2JqJ(l, dy/q)J(0, dy/q) Hq/L) 2 
-AqJJ{\, a^) H<l/L) + iqJJ(l,ay/q)J(0, dy/q) Hq/L) 

+4q 2 J(l, a^g) ln(g/L) - 4g 2 J(l, a^/g) J(0, a^g) ln(g/L)). (78) 
After applying a similar technique, it turns out that the 23 factors are 



(231) G „ 3 (2312)G . „ 4 (2313)G 

(232) G 3 (2323)G „ 4 (2324)G 

% = -P n s + 2p 77.^ , (79) 



where the 230 and 231 sub- factors are expressed as 



n4 230)G = 8Jp 3 y/qy/j - q(2 J(0, dy/q) + 2pJ(0, dy/q) ln(g/L) 

-pJ(0,aVg) 2 ln(g/i) - 7(0, a^) 2 - p!n(g/L) - 1), 

7^ 2312)G = Vj(2JqJ(l,dy/q)d + 2q 2 J(l,dy/q)dJ(0,dy/q) + 2Jy/^l(p) 2 

-2JqJ(l, a^/q)aJ{Q, dy/q) - 2q 2 J{l, ay/q)a + 2 Jy/qJ(0, ciy/q) 2 
~Zq 3 / 2 {{q/Ly) 3 - 3g 3 / 2 J(0,a^) 2 - 4J^J(0,a^) 
+6q 3 / 2 J(0,aV9)), 
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n (23i3)G = v /^j(2Jl n ( g /i) + 2JJ(0,a^) 2 ln(g/L) - 3gJ(0, a^/g) 2 ln(g/L) 
-3g ln(q/L) - 4 J J(0, a^/g) ln(g/£) - 4<? + 4 J + 4 J J(0, a^) 2 
-4<?J(0, aV?) 2 + 6gJ(0, a^g) ln(g/L) + 8<? J(0, a^g) 

-8JJ(0,oVg)). (80) 
While, the 232 factors come as per the following identifications 

n (2323)G = a f^J_ q ( q S/2 a J( 0j a ^f + J q 3 / 2 a J(0, dy/q) - Jq 3/2 aJ{0, Cly/q) 2 

-q 5/2 a.J(0, ay/q) + 2q 5/2 J(l, dy/qfa - 2q 3/2 JJ{l, dy/qfd 

-2q 2 J{l,a^)J{0,a^) - 2JqJ(l,dy/q) 

+2JqJ{l, dy/q)J(0, ay /q) + 2p 3 q 2 J(l, dy/q)), 
n (2324)G = a f^J— q ^/2j^^ a ^a\n{q/L) + q b l 2 aJ{0,a^f\n{q/L) 

-q 5/2 aJ{0,ay/q~) ln(q/L) - q 3/2 JJ(l,dy/q) 2 a\n(q/L) 

+Jq 3/2 aJ{0 1 dy/q) \n(q/L) - Jq 3/2 aJ(0, ay/q) 2 \n(q/L) 

+2JqJ(l, ayfq)J{0, dy/q) \n(q/L) - 2 JqJ(l, dy/q) ln(q/L) 

~2q 2 J(l, ay/q)J(0, dy/q) \n(q/L) + 4q 2 J(l, dy/q) 

+AqJJ{0, dy/q) J(l, dy/q) - Aq 2 J(l, dy/q) J(0, a^) 

+2g 2 J(l, dy/q) Hq/L) - 4q JJ(1, dy/q)). (81) 
We notice further that the factors of the ^(p) 3 -terms can be expressed as 



nf G = nf )G + (q/LYnf )G + {q/L) 2 Vnf )G + ((q/LfVnf )G ). (82) 
Following similar factorizations, we have the following terms 



(30)G 

s 


(300)G 
— ri s 


h4 30l ) G /- 


-4 302)G / 2 , 




(31)G 
S 


= 4 3l0)G - 


h4 3ll ) G /- 


-nf^f 2 , 




(32)G 
S 


(320)G 
— fig 


Vnf 1)G f- 


-4 322)G / 2 , 




(33)G 
S 


(330)G 

= n s 




-4 332)G / 2 , 


(83) 



which factor nicely. As per the expectation, we find that the 30-factors are 
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{300)G = 4y/q-y/J^7 l J(l-2J(0,ay/q-) + J(0,ay/q-) 2 ), 
n {301)G = 4aVjJ(l,ay/q-){q 2 -q 2 J(0,ay/q-)- Jq + Jqj(0,dy/q)), 

4 302 ) G = f 2 a 2 q 3 / 2 y/j—qJ(l,a^) 2 (J-q). (84) 



Similarly, the 31-factors are obtained as 



ii s 



= 12^F 1 y/J^pJ(l + J(0 1 ay/^) 2 -2J(0,a^)), 
(311)G = l2qpaVj(qJ(l,dy/q)-qJ(l,dy/q)J(0,dy/q) 



-JJ{l,dy/q) + JJ(1, ay/q)J(0, dy/q)), 
4 312)G = 3q 3 / 2 y/T^f 2 a 2 J(l,ay/q-) 2 p(J-q). (85) 



The 32-factors reduce to 
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Figure 16: The surface minor p G plotted as a function of the scale and angular momentum, viz., q, J, describing 
the nature of the stability of the qp-surface for general massive rotating quarkonia. 



f 0)G = l2^q^/j^p 2 J(l + J(0,a^q) 2 -2J(0,aVq)), 
= 12aqp 2 Vj(qJ(l,a^) - qJ(l,a^/q)J(0,a^q) 
+JJ(l, ay /q)J(0,a^q) - JJ(l,a^q)), 

(322)G 



(321)G 
l S 



» s - Zq A/2 ^J~qfa 2 J{l,a^fp 2 {J-q) 

Finally, it follows that the 33-factors are given by 



(86) 



(330) G 
n S 

(331) G 
S 

(332) G 
S 



I) 



II 



^^J^p 3 J(J(0,a^) 2 - 2JJ(0,a^) + 1), 
4aqp 3 V r JJ(l,a^/q)(q - qJ(0,a^/q) - J + JJ(0, a^/q)), 
+q 5/2 ^fJ~qfa 2 J{l,a^q) 2 p\j - q). 



(87) 



For the general quarkonia with J} fluctuating, the stability of the qp-surfacc requires that the principle 
minor p G remains positive on (M3, g). Correspondingly, this leads to the constraint that the {q,p, J} satisfy 



nf G + n^ G l{p)+nf G l{p) 2 +nf G l{pf < 0. 



(88) 



The three dimensional graphical views of the surface minor p G of generic quarkonia is depicted in the 
Fig. (|16|) . For the limiting small scale q — > + , the graphical view of the surface minor p G shows a large negative 
minima of the order of the height — 10 17 , which occurs in the regime of a large angular momentum J — > 5 — . 
When all the parameters arc allowed to fluctuate and b is an arbitrary real number, the exact intrinsic geometric 
characterization of the minor is explicitly obtained in the Eqn. (|59p , which forms the basis for the stability of 
the gp-surface. 

After considering the contributions of the fluctuations to the strongly coupled massive rotating quarkonia, 
we obtain the following general expressions for the determinant of the metric tensor 



P/(~1 + J"(0,a-v/g)) , ( o)G 

8bH(p) e exp {3l(p))q^ 2 J 3 / 2 ( J - qf/ 2 {ng 



n^ G l{p)+n^ G l(p) 2 +nf G l{pf 



(89) 
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where the ng°^ G -term factorizes as 



nf G = {{q/Lf ) 3 (8p 4 n<°*> G + 8p 5 n^ G + 8p 6 n^ G ) 



After simplification, we obtain the sub-factorizations 

n (04)G = n (040)G + /n (041) G) 
„(06)G = n (050)G + /n (051) G) 

n T )G = n T 0)G + /4 061)G 

A r (040)G (041)G (050)G (051)G (060)G (061)G-i , 

As per our computation, the {n g ,n g ,n g ,Ug ,n,g ,n,g j are expressed 



".</ 



U 9 



(04 0) G = V7yj^^ (2 j (0 , a ^)_ j(o, a ^) 2 -i), 



n (050)G = y^V J - 9\/g(4 J(0, aV?) ln(g/i) + 6.7(0, aV?) 

-2 ln( 9 /i) - 2 J(0, a^) 2 ln( g /i) - 3 - 3 J(0, a^) 2 ), 

n (051)G = 3g 3/2 + 3(7 3/2 J(0ia ^ ) 2 +2(? 3/2 J(0;a ^ ) 2 lll(g/L) 

+2g 3 / 2 hi(g/L) - 6g 3 / 2 J(0, a^) + 87^7(0, ajq) ln(q/L) 

+12Jy/qJ(0, aV?) - U^q\n(q/L) - 6JV«J(0, a^g) 2 

-4 J^J(0, a^) 2 ln(g/L) - 4g 3 / 2 J(0, aV?) Hl/L) - 6 J^, 



n (060)G = Jj^JJ—q^ j(o, ln( g /i) + 2 J(0, a^) ln((?/L) 2 
-2 ln(g/L) - 7(0, a^q) 2 \n(q/L) 2 - 2 J(0, a^) 2 ln(g/L) 
-ln(g/L) 2 ), 

n (061)G = g 3/2 ln((?/jL) 2 +(?3 /2 J(0;a ^ ) 2 ln((z/L) 2 + 2(7 3/2 ln((z/jL) 

+2g 3 / 2 J(0, a^) 2 ln(g/L) - 2 J^J(0, a^) 2 ln(q/L) 2 
-2g 3 / 2 J(0, aV?) M<zA) 2 - 4 9 3/2 J(0, a^/g) 
-4 JV? J(0, a^f \n{q/L) + 4 J^J(0, a ^ ln(g/i) 2 
-2 JV? ln(g/L) 2 + 8Jy/qJ(0, a^g) In(?/L). 

We find that the l{p) terms factorize as follows 

n^ G = {{q/LYfn^ G + ({q/Lffn^ G 
with the following sub-factorizations 

n (12)G = n ( 120 )G + /n (121)G j 

4 13)G = 4 130)G + M 131)G 

The factors of the 12-components are given by 
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4 12 ° )G = 4 1202) V + "< 1203) V + 4 1204 ) V, 

n (121)G = n (1212)G p 2 + n (1213)G p 3 + n (1214)G p 4 + n (1215)G p 5 (g7) 

with the following sub-factorizations 



n (1202)G = 2 VJy/j - q(a 2 q 3/2 J(0, a.y/q) 2 - a 2 q 3/2 J(0, ay/q) + 6J(l,a.y/q)aq 

-6J(1, a^)agJ(0, a^g") - 8^J(0, a^) 2 - 8^ + 16^J(0, a^)), 

n (i203)G = 4Vj\/J^~q(a 2 qV 2 J(0, a^f ln(q/L) - a 2 q 3 ' 2 J(0, a^g) ln(g/X) 

+6J(1, a,y/q)aq + 6J(1, ay/q)aqln(q/ L) — 6,7(1, ay/q)aqJ(0, ay/q) \n(q/L) 
-6?J(1, <Vg>J(0, aV?) - 8V9ln(g/i) - 8^9 J(0, a^) 2 MffA) 
-8^7(0, a^/g) 2 - 8V? + 16^J(0, a^g) + 16^(0, aVg) hi(g/L)), 

n (i204)G = 2 Vj v /7^(a 2 g 3/2 J(0, a^) 2 ln(g/L) 2 - a 2 g 3 / 2 J(0, a^) ln(g/L) 2 
-6J(l,a^q)aqJ(0, a^g) \n(q/L) 2 + 6J(l, a^)aq\n(q/ L) 2 + 8y/q~ 
-8y/qln(q/L) - 8J{l,ay/q)aqJ(0, a^g) ln(q/L) + 8y/qJ(0, a^) 2 
-8y/qln(q/L) 2 - 8y/qJ(0, a^) 2 \n(q/L) 2 - 8^J(0, a^) 2 \n(q/L) 
+8q hx{q/L)J{l, a^)a - 16^J(0, ay/q) + 16^qJ(0, ay/q) ln(g/L) 2 
+Wy/qJ(0,ay/q)ln(q/L)). (98) 

The 121-factors turn out to be 



n (i2i2)G = aq (2J q 1 / 2 J(l,ay/q) 2 a + 2aq 3/2 J(0 1 ay/q')-2aq 3/2 J(0,ay/q) 2 
+4gJ(l, av/5) J(0, a^) - 4J(1, a^)g + 4Jag 1/2 J(0, a^g) 2 
-4Jag 1/2 J(0, ay/q) - 8JJ(l,ay/q)J(0, ay/q) + 8JJ(1, ay/q) 
-q 3 / 2 J(l,ay/q-) 2 a), 

n (i2i3)G = 2(- g 5 / 2 J(l,a^) 2 a 2 ln(g/L) + 2a 2 g 5/2 J(0, a^g) ln(g/L) 
-2a 2 g 5/2 J(0, a^) 2 ln(g/i) + 2 Jg 3 / 2 J(l, a^)V ln(g/i) 
+4J(1, a^)ag 2 J(0, a^s) ln(g/L) + AJa 2 q 3/2 J(0, a^) 2 \n(q/L) 
-4g 3/2 J(0, ay/Fj) 2 - 47(1, a^g)ag 2 \n(q/L) - 6J(1, a^aq 2 
-4Ja 2 g 3/2 J(0, ay/q) hx{q/L) - Aq 3/2 + 6 J(l, a x /q)aq 2 J{Q 1 ay/q) 
+8Jy/q + 8g 3/2 J(0, Oy/q) - 8 J J(l, ay/q)aqJ{Q, ay/q) ln(g/Z) 
+8 J J(l, a^/g>g ln(g/Z) + 8Jy/qJ(0, a^) 2 + 12 J J(l, a^)ag 
-12JJ(l,aVg>gJ(0,Gh/g) - 16Jy/qJ(0,ay/q~)), 

n (i2i4)G = 2( _ 9 5/ 2 J(1) a ^ } 2 a 2 + 2 Jg 3 / 2 J(l, a^g) V ln(g/L) 2 

+2a 2 g 5/2 J(0, ay/q) \n(q/L) 2 - 2a 2 q 5/2 J{0, a^) 2 \n(q/L) 2 
+4J(1, ay/q)aq 2 J(0, a^s) ln(g/£) 2 - 4ag 2 J(l, ay^) ln(g/L) 2 
+4Ja 2 g 3/2 J(0, a^) 2 ln(g/L) 2 - AJa 2 q 3/2 J(0, a^?) ln(<?/ L ) 2 
+8 J(l, ay/q)aq 2 J(0, aVs) Mg/£) ~ 8 J ( 1 > a Vv) a( l 2 H<l/ L ) 
+8JaqJ(l, ay/q) ln(g/L) 2 - 8 JJ(1, a^)ag J(0, a^) ln(g/L) 2 
-16g 3/2 J(0, a^g) 2 ln(g/L) - 16 JJ(1, a^ag J(0, a^/g) ln(g/L) 
+16JJ(l,a x /9)«9 ln ('?/ i ) - 16g 3/2 ln(g/L) - 28g 3/2 J(0, a^/g) 2 
-28g 3/2 + 32 Jy/q]n(q/L) + 32 J^/g J(0, a^) 2 ln(g/L) + 56 J^g 
+32g 3 / 2 J(0, a^) ln(g/L) + 56g 3/2 J(0, a^?) + 56,7^(0, «\/g) 2 
-64J^J(0,a^)ln(g/L) - 112Jy/qJ(0,ay/q), 



35 



4 1215 > G = 4^\n(q/L)(-q\n( q /L)-qJ(0,a^) 2 ln(q/L) + 2JHq/L) 

-2q - 2g7(0, a^q-f + 2qJ(0, dy/q) Hq/L) + 2 7 7(0, a ^) 2 ln(g/L) 
-4 J J(0, a^/q) ln(g/L) + 4 7 7(0, a^) 2 + 4 J + 4g7(0, a^g) 
-877(0, aV5)). 

Similarly, the factors of 13-components are given as 

4 130)G = 4 1303) v + 4 1304) v + 4 1305) v, 

n ^ G = 4 1313 ) v + 4 1314) V + 4 1315) V + 4 1316) V- 

The 130-sub-factors turn out to be 

n U303)G = 2 <JjyJJ - q(a 2 q 3/2 J(0, a^qf - a 2 g 3/2 7(0, ay/q) + 67(1, d^ag 

-6ag7(l, dy/q)J(Q, ay/q) - 8y/q - 8^7(0, a^q) 2 + W^qJ(0, ay/q)), 

n (i304 )G = 4 V 7 v /ir^(oV / V(o,oV^ a Mg/i)-oV / V(o,oV5)iii(9/i) 

+6g J(l, a.y/q)a — 6ag7(l, ay/q)J(0, ay/q) ln(g/L) + 6 J(l, a^/q)aq ln(g/7) 
-6qJ(l,ay/q)aJ(p, ay/q) -8y/q- 8^qJ(0, a^f Hq/L) - 8y/q]n(q/L) 
-8y/qJ(0, a^q-) 2 + 16^(0, ay/q) Hq/L) + Wy/qJ(0, ay/q)), 
n (l305)G = 2 ^/j^T~q{a 2 q 3 / 2 7(0, ajq) 2 ln(g/L) 2 - a 2 g 3 / 2 7(0, ay/q) ln(g/L) 2 

+6 J(l, aV?)ag hi(g/L) 2 - 6 J(l, aVg)ag7(0, a^/g) hi(g/L) 2 - &y/qhi(q/L) 
+8^,1(0, a^) 2 - 8^7(0, a^f \n(q/L) 2 + 8^- 8^\n(q/L) 2 
-8qJ(0, ay/q) ln(g/7)7(l, a^)a + 8q\n(q/L)J(l,a,y/q)a 
-8y/qJ(0, a^) 2 ln(g/L) + 16^7(0, a^g) ln(g/L) 2 
+ 16^7(0, a^g) Hl/L) - 16^7(0, 0^/5))- 

The 131-factors take the following explicit expressions 

n (l3i3)G = a q{-q 3/2 J{l,a^) 2 a + 2Jq 1/2 J{l,a^) 2 a + 2aq 3/2 J{Q,ay/q) 
-2ag 3 / 2 7(0, a^q) 2 - 4J(1, a^)g + 4J(l,ay/q)qJ(0, ay/q) 
-47ag 1/2 7(0, ay/q) + 47ag 1/2 7(0, a^q) 2 + 8 JJ(1, a^g) 
-877(1, aV5)7(0, 0^5)), 

4 1314 ) G = 2(-g 5 / 2 7(l, a ^) V i n ( g / L ) _ 2 g 3 / 2 - 2g 3 / 2 7(0, a^) 2 

+2a 2 g 5/2 7(0, a^g) ln(g/Z) - 2a 2 g 5/2 7(0, ay/q) 2 Hq/L) 
+27g 3/2 7(l, ay/qfa 2 ln(g/7) + 47^9 + 47^7(0, a^g") 2 
+47a 2 g 3 / 2 7(0, aV?) 2 ln(g/7) + 4ag 2 7(1, a^) 7(0, aV?) M<?/£) 
-47(1, a^/q)aq 2 ln(g/7) - 47a 2 g 3/2 7(0, a^g) ln(g/7) 
+4g 3 / 2 7(0, oV5) - 67(1, <Vg>g 2 + 67(1, a% /g>g 2 7(0, a^q) 
-87^7(0,0^/9) + 877(1, aV5)agln(g/L) + 1277(1, a^ag 
-87ag7(l,a 1 /g)7(0,Gh/g) ln(g/L) - 1277(1, a 1 /g)og7(0, Gh/g)), 
4 1315 > G = -g 5 / 2 7(l,aV5) 2 « 2 ln(g/i) 2 + 2a 2 g 5 / 2 7(0,a^)ln(g/L) 2 

-2a 2 g 5 / 2 7(0, a^/g) 2 ln(g/L) 2 ) + 27g 3/2 7(l, a^) 2 a 2 ln(g/L) 2 
+47(1, aVg>g 2 7(0, ay/q) ln(g/L) 2 - 47(1, <Vg>g 2 ln(g/L) 2 
+47a 2 g 3/2 7(0, a^) 2 Hq/L) 2 - 47a 2 g 3 / 2 7(0, ay/q) hi(g/i) 2 
+877(1, ciy/q)aq ln(g/7) 2 - 877(1, CLy/q)aqJ(0, ay/q) ln(g/7) 2 
+87(1, ay/q)aq 2 J(0, ay/q) ln(g/i) - 87(1, ay/q)aq 2 ln(g/7) 
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~8q 3/2 J{0,a^) 2 ]n(q/L) - 8q 3/2 \n(q/L) + 16Jy/qln(q/L) 
+16q 3/2 J(0,a^) ]n(q/L) + 16 JJ(1, ay/q)aq \n{q/L) - 16q 3/2 
-16JJ(1, ay/q)aqJ{0, a^) ]n(q/L) + 16Jy/qJ(0, ay/q) 2 \n(q/L) 
-I6q 3/2 J(0, a^) 2 + 32Jy/qJ(0, a^) 2 + 32g 3 / 2 J(0, ay/q) 
+32 - 32 J^J(0, ay/q) \n{q/L) - 64J^J(0, ay/q), 
4 1316 ) G = -32J^ln(g/i). (102) 

It follows that the l(p) 2 terms factorize as 

nf G = (q/Lyn^ G + ((q/L)n\^ G + ((q/Lrfn^ G , (103) 
where the sub-factorizations take the following forms 



n (21)G = n (210)G + /n (211)G ) 
n (22)G = n (220)G + /n (221)G 



h G = 4 23o)g +/4 : 

Explicitly, we find that the individual factors appear as 



n (23)G = n (230)G + /n (231)G^ (1Q4) 



(210)G = n (2101)G J3 + n (2102)G p 2 + n (2103)G p3) 



n (211)G = n (2111)G p + n (2112)G p 2 + rl (2113)G p 3 (1Q5) 

The 210-factors are given by 



4VJy/j- q(a 2 q 3/2 J(0, ay/q) - q 3/2 J(l, ay/q) 2 a 2 - a 2 q 3/2 J(Q, ay/q) 2 
-2,7(1, a^q)aq + 2J(1, ay/q)aq J(0, a^) +4^9 + 4^J(0, ay/q) 2 
-8y/qJ(0,ay/q))), 

4\[j\/J - q{a 2 q 3/2 J(0, ay/q) Mq/L) - a 2 q 3/2 J(0, ay/q) 2 Mq/L) 
+6 J(l, ay/q)aqj(0, ay/q) ln(q/L) + Qp 2 J(l,ay/q)aqJ(0, ay/q) 
—6J(l,ay/q)aq — 6J(l,a^/q)aq\n(q/L) + 8^/q\n(q/L) 
+8y/qJ(0, a^) 2 Mq/L) + 10V? + 10^(0, ay/q) 2 
-16V?J(0, ay/q) Hq/L) - 20^J(0, ay/q)), 

2VJ^J - q{q 3 ' 2 J{l,ay/q) 2 a 2 Hq/L) 2 + 4/(1, a^aq J(0, ay/q) \n(q/L) 2 

+4qJ(0, ay/q) ln(g/L) J(l, a^)« + V5 J (0. aV?) 2 + Wn^/L) 2 
-4gln(g/L) J(l, aVg)a - 4.7(1, a^)ag ln(g/£) 2 - 8y/qJ(Q, ay/q) 
+4^9 + 4^(0, aV?) 2 ln(g/L) 2 + 8Vgln(?/L) + 8^J(0, a^) 2 \n(q/L) 
-8^J(0, aV?) Hq/Lf - 16^(0, a^g) In(g/L)). (106) 

While, the 211-factors are given as 

„(ani)G = 2a(a g 5 / 2 J(0,a^) 2 -a g 5/2 J(0,aV7i) + g 5 / 2 J(l, a 77i) 2 a 

+2 J(l, a^)g 2 - 2 J(l, a^)q 2 J(0, a^) - 2 Jg 3/2 J(l, aV?) 2 « 
-2 J a(? 3/2 J(0, a^Tj) 2 + 2Jpaq 3/2 J(0, ay/q) - 4JJ(1, a^g 

+4Jj(i,aVg)gJ(o,aVg)), 

n (2ii2)G = 2a( g 5 / 2 J(l,a^) 2 aln( (? /L) + 2 a q 5 / 2 J(0, a ^) 2 ln( (Z /i) 



(2101)G 



n (2W2)G 



(2103)G 
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-2Jq 3/2 J(l, av ^) 2 a\n(q/L) - 2aq 5/2 7(0, ay/q) ln(g/L) 
+4 Jag 3 / 2 J(0, a^g) In(g/7) + 4J(1, a^)g 2 ln(g/Z) 
-47(1, a^g^O, a^g") ln(g/L) - 47ag 3/2 J(0, a^g) 2 ln(g/L) 
+6g 2 J(l,a^) - 6g 2 7(0,a^)7(l,a^) - 12JqJ(l, ay/q) 
+8JJ(l,ay/q)qJ(0, ay/q) Hq/L) - 877(1, ay/q)q ln(g/i) 
+12JqJ(0,ay/q)J(l,ay/q)), 
nf 13 ) G = aq(aq 3 / 2 J(0 1 a^) 2 Hq/L) 2 + 2qJ(l 1 a^)Hq/L) 2 

-2J(1, aV?)gJ(0, aV?) ln(g/£) 2 - ag 3/2 7(0, a^g) Hq/L) 2 
-4g7(0, a^g") ln(g/L)7(l, ay/q) + 4gln(g/Z)7(l, a^g) 
+8 JJ(0, a^/q-) ln(g/L)7(l, ay/q) - 87ln(g/Z)7(l, ay/q) 
-4JJ(1, aVs) Hl/L) 2 + AJJ(l,ay/q)J(0, ay/q) ln(g/L) 2 

+27ag 1/2 7(0, ay/q) ln(g/L) 2 - 27ag 1 / 2 7(0, ay/q) 2 Hq/L) 2 ). (107) 
The 22-factors can be cascaded into the following form 



4 220)G = 4 2202) V + 4 2203) V + 4 2204) V, 



B 

with the individual sub-factors 



n (221)G = n (2212) Gp 2 +n (2213) Gp3+rl (2214) Gp 4 (10g) 



n 



(2202)G 



(2203)G 



n 



(2204)G 



2^/1 \J J - q(U(l,ay/q)aqJ(0, ay/q) - Aa 2 q 3/2 7(0, a^?) 2 
-4J(1, <Vg>g + 4a 2 g 3/2 7(0, a^/g) - 5g 3/2 7(1, ay/q) 2 a 2 
+12^+12^7(0, a^g") 2 - 24^7(0, ay/q)), 

4V7 v / 7-g(-g 3/2 7(l, a^f) 2 a 2 ln(g/L) + 2a 2 g 3/2 7(0, a^/Tf) ln(g/L) 
-2a 2 g 3/2 7(0, ay/q) 2 ln(g/i) + 87(1, ay/q)aqJ{Q, ay/q) ln(g/£) 
-87(1, ay/q)aq ln(g/L) + 107(1, ay/q)aqJ(Q, ay/q) 
-10qJ(l,ay/q)a + 12Vg7(0, a^Tf) 2 ln(g/L) + 12Vgln(g/L) 
+ 16V? + 16Vg7(0, aV?) 2 - 24^7(0, a^) ln(g/L) - 32^7(0, ay/q)), 
2V7v/7-g(g 3/2 7(l, a^)V ln(g/L) 2 + 4g7(0, a^g) ln(g/L)7(l, a^)a 
+4^g - 47(1, ax /g)agln(g/i) 2 + 47(1, ay/q)aqJ(0, a^g) ln(g/L) 2 
-4gln(g/L)7(l, + 4Vg7(0, a^f ln(g/L) 2 + 4,/g ln(g/L) 2 

8^7(0, a^)ln(g/L) 2 
-16Vg7(0,aVg)ln(g/£)). 



+4 v /7-g7(0,a^) 2 - 8^7(0, a^g) - 
+8Vgln(g/i) + 8^7(0, a^f Hq/L) 

The 221-factors are given by the following expressions 



(109) 



n (2212)G = a ( 3a? 5/2j(0, _ 3aq 5/2j^ + ^5/2 a ^f a 

+67ag 3 / 2 7(0, o^/g) - 67(1, a^)g 2 7(0, ay/q) + 67(1, a^g 2 
-67ag 3/2 7(0, a^) 2 - 87g 3/2 7(l, ay/q) 2 a - 127g7(l, a^g) 

+1277(1, aVg)g7(0,aVg)), 

n (22i3)G = 2ag(2g 3/2 7(l,a x /7i) 2 aln(g/L) + 3ag 3/2 7(0,a x /7i) 2 ln(g/ J L) 
-3ag 3/2 7(0,a^) ln(g/L) - 47g 1/2 7(l, a^) 2 aln(g/L) 
+67(1, ay/q)qHl/L) - 67(1, a^)g7(0, ay/q) ln(g/i) 
+67ag 1/2 7(0, a^/g) ln(g/7) - 67ag 1/2 7(0, a^/g) 2 ln(g/L) 
-10g7(0, aVg)7(l, ay/q) + 1277(1, ay/q)J(0, ay/q) ln(g/Z) 
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+10?J(1, oVg) - 12JJ(l,Oy/q) Hq/L) - 20 JJ(1, ch/q) 
+2DJJ{0,a y /q)J(l,( h /q)), 
4 2214)G - aq(aq^ 2 J(0, a^)' ^A) 2 - a<? 3/2 J(0, aV?) W^) 2 

-2 J(l, ay ^)qJ(0, Oy/q) ln(g/L) 2 + 2Jaq 1/2 J(0, a^/q) \n(q/L) 2 
-2Jaq 1 ' 2 J(Q, a^q) 2 \n(q/L) 2 + 2qJ(l, a^) \n(q/L) 2 
+Aq ln(q/L)J(l, ajq) - Aq J(0, a^/q) ln(q/L)J(l, a^/q) 
-AJJ(l,ay/q) Hq/L) 2 + AJJ(l,a^q)J(0, a^/q) ln(g/L) 2 

+8JJ{0,a^)ln(q/L)J(l,a,/q) - &Jhi(q/L)J(l,a,/q)). (110) 
Similarly, the 23-factors turn out to be 



n ( 23 °) G = 4 2303 'V + 4 2304) V, 

n (231)G = n (2313)G p 3 + n (2314) G;A (ln) 



n (2303)G = 2\f~J^fq\J J - q{2a 2 qJ(0 1 a^/q) - 2a 2 qJ(0, a^/q) 2 - 3qJ(l, a^fqfa 2 
+4J(0, a^) 2 + A - 8 J(0, a^/q)), 

n (2304)G = 4 V7 v ^7^( a 2 (7 3 / 2 J(0, ay /q) ln(q/L) - a 2 qV 2 J(0, a^q-) 2 ln(q/L) 
-q 3/2 J(l, a^/qfa 2 hi(q/L) + 2aqJ(l, a,y/q)J(0, a^/q) Hq/L) 
-2aqJ(l,ay/q) ln(q/L) + A^Hq/L) + AJ(l,a,y/q)aqJ(0, a^q) 
+4^J(0, a^q) 2 HQ/L) - 47(1, a^)aq + 6^ + 6^J(0, a^q) 2 
-8^J(0, a^) H^/L) - 12y/q J(0, Oy/q)), (112) 



4 2313 ) G = aq(aq 3 / 2 J(0 1 a^) 2 -aq^ 2 J(0,a^) + 2 q ^ 2 J(l,a^) 2 a 
+2qJ(l,<h/q) - 2 J(l, a^)qj(0, a^q) + 2Jaq 1/2 J(0, a^) 
-2Jaq 1/2 J(0, ay/q) 2 - AJq 1/2 J(l,a^/q) 2 a - 4JJ(1, a^) 
+4JJ(l,a^)J(0,a^)), 

n ^23i4)G = 2a g (a 9 3/2 J(0,a^) 2 ln( g /L) + ( 7 3/2 J(l,a^) 2 aln( g /L) 
-aq 3/2 J(0,ay/q) ln(g/L) - 2Jg 1/2 J(l,aVg) 2 alu(<?/£) 
-2 Jag 1/2 J(0, a^g) 2 ln(g/L) - 2 J(l, a^Jq)qJ{Q, a^/q) ln(g/L) 
+2gJ(l, a^s) ln(g/i) + 2 Jag 1 / 2 J(0, a^) ln(g/L) + 4gJ(l, a^s) 
-4?J(0, a^/g) J(l, aV?) + 4JJ(1, Oy/q)J(0, a^) ln(g/X) 

-4JJ(1, o^/g) ln(g/L) + 8 JJ(0, ay /q)J{l,ch/q) - 8 JJ(1, a^g)). (113) 
Finally, the l(p) 3 terms are given by 

n (3)G = n (30)G + {q/L )P n mG + {{q/L)P) 2 n (32)G + {{q/ L) Vf n (33)G ^ (m) 

where the factors are expressed as 



(30)G 

'9 


= n (300)G _ 


h/4 3oi ) G , 


(31)G 

'9 


= nf°) G - 


h/»< 311)G > 


(32)G 

'g 


= 4 32 °) G - 


h/4 32i > G 


(33)G 
9 


= 4 33 °) G - 


h/4 33i ) G . 



(115) 
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Figure 17: The dimensional behavior of the determinant of the metric tensor \\g\\ plotted as a function of the 
scale and angular momentum, viz., q, J, describing the nature of the stability of the general massive rotating 
quarkonia. 



The 30-factors are given by the following expressions 



„(300)G = 2aqVJy/J - q{aq 1 l 2 J{Q 1 a^) 2 - aq 1 ' 2 J(0, ajq) 



+27(1, ay/q) J(0, a^/q) + 2q 1 ' 2 J{l, a^) 2 a - 2 J(l, a^q)), 



n (3oi )G = a v /2 J(i,aVg) 2 (2J-g). (116) 

The 31-factors take the following forms 

nf w)G = Qapq^fjyjj - q(aq 1/2 J(0, a^/q) 2 - aq 1/2 J(0, a^/q) 

-2J(1, a^q) + 2J(0, a^q)J(l,a^q) + 2q 1 / 2 J(l, a^) 2 a), 
n (311)G = ^ q ^J^ a ^qf p{ 2J-q). (117) 

Similarly, the 32-factors turn out to be 

n (320)G = QaqpiyTJ^j _ q (aq 1/2 J(0, a^f - aq 1 ' 2 J(0, a^/q) 

-2 J(l, aV?) + 2q^ 2 J(l,a x /qfa + 2.7(0, a^g) J(l, <V«)), 
n (32i)G = 3a Vq 3 / 2 J(l, a ^) 2 (2J- g ). (118) 

Finally, the 33-factors can be obtained as follows 



n 



f°) G = 2agp 3 ^ v ^r^( ag 1 / 2 J(0,a^) 2 - ag 1 /2 J(0;a ^ ) 



+27(1, a^) J(0, ay/q) - 2 J(l, a^) + 2q 1 / 2 J(l, a^) 2 a), 



n 



< 331 ) G = a yj(l, a V9)V(l,<V7i) 2 (2J- g ). (H9) 



From the three parameter quarkonia, the Fig. (|17[) offers the general three dimensional behavior of the 
determinant of the metric tensor. From the viewpoint of the fixed p quarkonia in the Regge regime as previously 
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described, we notice in the present case that the graph of the determinant of the metric tensor shows the similar 
nature of the Gaussian fluctuations of the strong QCD coupling. Notice that the underlying quarkonium 
configuration may become unstable in specific regions of the parameters. For instance, for small values of the 
momentum scale and angular momentum, viz q, J — > + , we observe that the determinant of the metric tensor 
picks up a large amplitude of order 10 21 . However, for intermediate values of q and angular momentum q, J, the 
Fig. fTT)) shows that the underlying configuration is well-stable and so is the basis for the three dimension intrinsic 
manifold (M 3 ,g). From the observation of the Fig. (JTTJ) , we predict that the regions of the thermodynamic 
stability are present in the interval q, J € (1,4). Globally, the stability of (A/3, g) constraint the principle 
minors {gu,p G , \\g\\} to remain positive. Specifically, for the same sign of {b,p, /}, the volume stability of the 
(A/3, g) imposes the following constraint 

nf G + n^ G l(p)+n^ G l(pY + nf G l{pf > 0, J(0, arf) < 1, 

< 0, J(0,Oy/q) < 1. (120) 

Importantly, it is worth mentioning that both the limiting configurations with J = q and J(0, a^Jq) = 1 arc 
abided from the thermodynamic stability constraints. To summarize the phases of generic quarkonia, the exact 
formula for the scalar curvature may analogously be deduced as the one we have offered for the q J-planc. With 
some modification of the pre-factors, we find that the summation over l(j>) naturally arises with the B n as the 
polynomials in p, whose coefficients can be expressed as the functions of the Bloch-Nordsieck logarithmic factor 
f(q,J). In the present case, the denominator of the scalar curvature involves the third and second power of 
the factors in the numerator of the determinant of the metric tensor. Interestingly, the quantitative property 
of the scalar curvature follows, except the fact that its explicit presentation is a bit prolonged. Up to a phase 
of QCD, we observe that the global properties of the three parameter quarkonia remain the same as we have 
exactly indicated for the Q J- and QM-planes. 

5 Discussion and Conclusion 

We have studied the thermodynamic geometric properties of non-abelian quarkonium bound states arising 
from the consideration of the one-loop gluon confining QCD potentials. Following the Polyakov argument, the 
thermodynamic intrinsic geometric nature of the QCD coupling is analyzed for the case of the Coulambic and 
rising confining regimes. Without any approximation, the effect of Regge and Bloch-Nordsieck rotations is 
explored form the stability perspective of the strongly coupled quarkonia with and without the inclusion of the 
effects coming from the mass. The fact that once the quark matter being formed should be stable, motivated us 
to determine the values of the parameters of strongly coupled quarkonia. We have obtained the allowed range of 
parameters from the stability conditions of the thermodynamic geometry. We have initiated our analysis with 
the consideration of massless non-rotating quarkonia. In this case, we have derived that the thermodynamic 
stability constraints impose the requirements ga > and \\g\\ > 0. Given the momentum scale q and index p 
of the configuration, our analysis determines an intrinsically stable inter dependence of the parameter from one 
and the other. The Regge rotating quarkonia is geometrically trivial because it only requires the positive sign 
of the rotation term. With and without the Regge rotation, the domain of the index p(q) can be offered as 

V m= v := {p I g lt > 0, \\g\\ > 0}. (121) 

In the case of massive rotating quarkonia, there is a set of thermodynamic stability constraints on the 
parameters which arises from the fluctuations of the q-angular momentum plane, q-mass plane and the entire 
q-angular momentum-index manifold, when all the parameters of the configuration are allowed to fluctuate. It 
turns out that one can fix all the parameters of the configuration following from one constraint to the other. 
One such procedure is to respectively consider the fluctuations on the QJ-plane, then on the QAZ-planc and 
finally on the QpJ-manifold. In the case of the Q J-plane, it follows that we determine the thcrmodynamically 
stable Q 2 (J) directly from the fluctuations of the Bloch-Nordsieck rotation. By taking an account of the 
thermodynamic stability constraint of the Q 2 (J), we can subsequently determine the thermodynamically stable 
domain of the mass M(Q) on the QM-plane. Finally, the thermodynamically stable domain of the index p can 
be determined as the function of Q and J on the Qp J-manifold. Thus, the consideration of the thermodynamic 
fluctuation theory offers the determination of the stable domains of the QCD strong coupling. This follows from 
the fact that the stable domain of q and the Bloch-Nordsieck rotation can be determined from one another, 
and thus the index can solely be expressed as a function of one the either two. This shows the power of the 
thermodynamic stability condition that determines all possible stable phases of the strongly coupled quarkonia. 
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Based on the analysis of the present paper, the thermodynamically stable index of the Bloch-Nordsieck rotating 
massive quarkonia is constrained by the following set 

Vrr^o := {p I n% > 0, ng > 0, nf > 0}. (122) 

To be specific, we have introduced the intrinsic geometric notion to QCD thermodynamics, and thereby it 
has been applied to study the behavior of constant and variable index strong QCD coupling in several simple 
thermodynamic considerations. Our explicit computations of the local and global thermodynamic corrections 
provide a prominent realization of the limiting equilibrium thermodynamic configuration of the strongly coupled 
QCD. This offers an intrinsic geometric exercise of the quarkonia, which has been described by the parame- 
ters of the strong coupling potential. The scalar curvature of the underlying thermodynamic metric tensor 
takes an exact form, which signifies the correlation volume of the underlying quarkonia. With and without the 
variable index strong QCD coupling, our consideration of the underlying thermodynamic fluctuations indicates 
interesting relationships between the geometrical concepts of quarkonium thermodynamics and phases of the 
massless, Regge rotating and Bloch-Nordsieck rotating quarkonia. For the case of the Bloch-Nordsieck rotating 
quarkonia, we have implemented the above constraints one by one, and thereby we conclude that the thermo- 
dynamically stable domain of the parameters is highly constrained on the (M^^g). For the value of p = 5/6 
which generically corresponds to the Regge rotating quarkonia, we observe that all the stability constraints are 
satisfied, whenever g,J £ (1,4). Thermodynamically, this determines all possible intermediate stable phases 
of the strongly coupled quarkonia. As the gluonic interactions become softer and softer, we find in the limit 
of Bloch-Nordsieck resummation that the underlying Sudhakov form factor offers the thermodynamic stability 
properties of the w, K and D s particles. 

Mathematically, we have computed the intrinsic geometric properties of the Besscl function of first kind con- 
voluted with two logarithmic functions of the respective weights (0, 2p). Explicitly, we determine the domain of 
the convolution of the zeroth order Besscl function J(0, a^/q), where the configuration finds globally stable fixed 
points and remains globally regular on the intrinsic three manifold (Ms, g). Notice that the exact determination 
of the global stability set V m jto is consistent with the local stability constraints on (M 3 ,g). Geometrically, we 
find that the global stability requirement of (M 3 ,g) avoids the limit J(0,a,y/q) — > 1. From the perspective of 
the intrinsic Riemannian geometry, it would be worth to find the globally stable and globally regular domains 
of the convoluted Bessel functions with finitely many elementary functions of distinct weights. This study 
herewith anticipates further examination of the globally stable and regular mapping class convolutions from the 
perspective of intrinsic Riemannian geometry. 

Physically, it would be interesting to analyze whether this kind of approach can be pushed further to explore 
intrinsic geometric properties of strong coupling QCD thermodynamic correlations. Specifically, we would like 
to understand the modifications towards the elements of the set P m ^o an d the associated Ricci scalar curvature, 
which we have explicitly offered for the Q J-plane. This geometric initiation would take a better shape once we 
take finite size effects of higher order perturbative QCD into account, and would reach its final goal with the 
complete account of non-perturbative QCD. The geometric feature of these explorations is left for the future. 
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